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ABSTRACT 


Behavior of continuous composite beams is investigated 
theoretically and experimentally. Analyses for local flange buckling 
and lateral buckling in negative moment regions are proposed based on 
torsional buckling theory and thin-walled beam theory, respectively. 
Deformations, including the effects of shear and slip, are studied in 
the elastic and inelastic regions. Procedures based on Newmark's in- 
tegration, finite difference and virtual work are employed for the eva- 
luation of moment-curvature relationships, beam deflection due to bend- 
ing and shear deformation, respectively. Three two-span continuous 
composite beams were tested. The test beams varied in terms of steel 
section size and amount of longitudinal slab reinforcement in the nega- 
tive moment region. Behavioral studies are conducted for failure loads 
and failure modes for a series of beams in which steel section size, 
concrete slab thickness and amount of longitudinal slab reinforcement 
are varied. Analyses satisfactorily predict ultimate loads, failure 
modes and behavior throughout the loading range. Based on experimental 
and analytical results design requirements are proposed for ultimate 
strength design of continuous composite beams. 

It is concluded that the maximum flange width-thickness 
ratio for the steel section should be less than 54/ yo if the amount 
of longitudinal reinforcement in a negative moment region is greater 


than the web area of the steel section. Strength of shear connectors 
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in a negative moment region is comparable with the strength of connectors 
in a positive moment region if sufficient development length of the 
longitudinal slab reinforcement is provided into the positive moment 


region. 
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CHAPTER I 
INTRODUCTION 


1.1 Introductory Remarks 


Presently in North America composite design is based to a 
large extent on elastic theory (1) (2) | However elastic design is not 
entirely routine in the case of continuous beams, because the beam stiff- 
ness is different in positive and negative moment regions. An ultimate 
design approach would eliminate this complication. Design based on 
ultimate strength may tend to produce a reduction in steel section size 
and in total depth. This reduction is related to the shape factor of the 
composite section which may reach a value of 2.0 in a positive moment 
region and 1.6 in a negative moment region. These values are signifi- 
cantly larger than those for wide flange steel sections. 

Extensive research on ultimate strength of simple span beams 
has resulted in the first requirements for limit design of simply-supported 
composite beams published in Britain in 1963 as Part I of CP 117°), How- 
ever further research on ultimate strength of continuous composite beams 
is still required, especially concerning conditions in a negative moment 
region. Research at the University of Alberta has included investigations 


by Piepgrass ‘*?, Teva Deven? into the behavior of composite 


beams under negative bending. 
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1.2 Previous Research on Continuous Composite Beams 


Research on ultimate strength of continuous composite beams 
has been conducted since 1960 at several universities including the Uni- 
versity of Cambridge, Imperial College, Lehigh University, University of 
Canterbury (New Zealand), University of Warwick, and University of Alberta. 
A brief summary of main conclusions of this research is presented in this 
section. 


(7) 


Culver et a tested a two-span continuous beam with a 
small amount of longitudinal slab reinforcement to establish the feasi- 
bility of designing continuous composite beams on the bases of ultimate 
strength. The results of the investigation concluded that only the steel 
section is effective in the negative moment region, and an expansion 
joint in the slab should be provided in the negative moment region to 
permit large rotation. 

Barnard (8) (9) tested a series of four continuous three-span 
beams in order to study the behavior of the most highly strained critical 
sections and the effect of their combined action on the behavior of the 
beam aS a whole. Two beams failed by crushing of concrete in the posi- 
tive moment region and the other two failed due to lateral buckling of the 
steel section in the negative moment region. The results indicated that 
continuous composite beams can be designed by the conventional plastic 
hinge method, if the positive moment hinge is the last to form; that the 
presence of slip strain does not significantly change the moment-curvature 
relationships from those obtained by assuming no slip strain; and that 
longitudinal slab reinforcement strengthens the negative moment regions, 


provided that buckling and shear failure can be avoided. 
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(10) conducted tests on fifteen 


Johnson, Van Dalen and Kemp 
Simply-supported beams in negative bending and six two-span continuous 
beams. The continuous beams were loaded in such a manner that they were 
affected by transverse bending in the slab. The slabs were reinforced 
with various amounts of longitudinal reinforcement. The results indi- 
cated that complete interaction occurred in negative moment regions 
where the concrete slab was cracked due to longitudinal tension. The 
actual ultimate loads for the continuous beams exceeded theoretical values 
by 11 to 31 percent. 


(11) tested to failure four continuous 


Daniel and Fisher 
beams which had been previously fatigue tested. The beams contained 
longitudinal reinforcement in the negative moment region. In each beam 
crushing of the concrete at the load points occurred at maximum load 
except in the beam with the largest amount of reinforcement where local 
buckling also occurred at ultimate load. Furthermore the longitudinal 
slab reinforcement in the negative moment regions attained its yield 
Stress at ultimate load. It was concluded that plastic analysis ade- 
quately predicted ultimate loads. 

park (12) tested four two-span continuous beams with various 
amounts of longitudinal slab reinforcement in the negative moment region. 
It was concluded that continuous beams with longitudinal slab reinforce- 
ment in the negative moment regions can be satisfactorily analyzed by 
simple plastic theory. Very little redistribution of bending moments 
was required in the test beams to develop an ultimate condition, since 
there was little difference between the distribution of elastic and ul- 


timate bending moments in spite of the large variations in the longi- 
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tudinal reinforcement content. Shear connectors, designed on the basis 
of ultimate capacity, performed adequately. 

Extensive research on the behavior of stud shear connectors 
in a negative moment region was conducted by Van Dalen et a1 (10) (13) (14) 
As a result, it was recommended that the design load for connectors in 
a negative moment region should be taken as 80 percent of the design 
values in a positive moment region. 


In 1970 C1imenhaga (1°) (16) 


investigated theoretical and ex- 
perimental post-buckling behavior of composite beams in negative moment 
regions. His proposed approximate method, based on the upper-bound plastic 
limit theorem, gave a reasonable prediction of the moment-rotation charac- 
teristics for I-shape sections after local buckling. 


(17) (18) (19) tested beams connected to columns 


Daniel et al 
simulating part of a multi-story frame. It was concluded that the ulti- 
mate strength of a composite section adjacent to the column face and su- 
jected to positive bending can be conservatively based on the strength of 
the steel section plus a portion of the concrete slab whose width equals 
the column width. Results indicated that ultimate strength of a composite 
section in negative bending can be based on the steel section plus the 
longitudinal slab reinforcement. Test beams exhibited sufficient rota- 


tion capacity in the vicinity of the joint so as to develop a mechanism 


at ultimate load. 


1.3 Previous Investigations at the University of Alberta 


Research on continuous composite beams at the University of 


Alberta was initiated in 1965 by Ferrigno with behavioral studies of 
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composite beams under positive moment. The slab dimensions used in 
Ferrier's test beams were chosen so that in one beam the neutral axis 
was in the slab whereas in the other the neutral axis was below the 
Slab. It was found the beams deflected more at working loads than 
predicted by usual analysis. This difference was likely due to resi- 
dual stresses in the steel section; however residual stress did not 
influence the ultimate strength. Ultimate concrete slab strains in ex- 
cess Of 0.005 were developed. Slip between the concrete and steel sec- 
tion did not prevent the beams from developing their theoretical ultimate 
moment values, provided that the total strength of the shear connectors 
supplied the required horizontal shear resistance. 


(4) 


The second investigation, by Piepgrass in 1968, was con- 
cerned with the behavior of composite beams in an isolated negative 
moment region. The main variables in this investigation were the amount 
of longitudinal reinforcement and slab width. Compression flanges of 
some steel sections were reinforced with cover plates in order to pre- 
vent premature local buckling. Beams without cover plates failed in 
local buckling in the flange and web whereas beams with cover plates 
failed in lateral buckling. Ultimate moments ranged from 98 to 107 
percent of the theoretical plastic moment values. The average stress 

in the longitudinal slab reinforcement at ultimate moment was less than 
the yield stress, which implied little rotation capacity. It was there- 
fore suggested that the use of steel sections defined as compact under 
the provisions of CSA Standard S16 may not necessarily guarantee adequate 


rotation capacity for the formation of a plastic hinge in a composite 


beam in negative bending. 
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Extensive tests of composite beams under negative moment 
were conducted by Parco! in 1969 and taal in 1970. The main 
variables in these tests were size of steel section, amount of longi- 
tudinal slab reinforcement and slab width. All beams failed in local 
buckling of the compression flange. It was concluded that significant 
increases in the negative moment capacity of composite beams can be 
achieved by the addition of longitudinal reinforcement, but these in- 
creases are not directly proportional to theoretical simple plastic 
moment values. For a given steel section, an increase in amount of longi- 
tudinal slab reinforcement results in a significant reduction in the ro- 
tation capacity. Therefore, to ensure the formation of a mechanism in 
beams where the negative hinge is first to form, it may be necessary to 


limit the amount of longitudinal slab reinforcement. 


1.4 Object and Scope 


This dissertation consists of both theoretical and experi- 
mental studies of the behavior of continuous composite beams. The ex- 
perimental study provides further information on deformations, moment 
redistributions and failure modes for continuous beams. Since relatively 
little analytical investigation has been conducted into local and lateral 
buckling in negative bending, an objective of this investigation is a 
theoretical determination of ultimate strength and rotation capacity in 
a negative moment region. A further objective is the prediction of 
ultimate load and failure modes for continuous composite beams. The 
study also includes a theoretical analysis for the deformation of con- 
tinuous composite beams. Finally design criteria for continuous com- 


posite beams are proposed. 
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CHAPTER II 


STRENGTH OF COMPOSITE BEAMS 


2.1 Introduction 


Characteristics of composite sections are different in 
positive and negative moment regions. In a positive moment region the 
concrete slab or a portion thereof is in compression and contributes 
Significantly to the moment resistance of the section as in Figure 2.1(a). 
Behavior under positive moment therefore is related to characteristics 
of the steel section, concrete slab and shear connectors. 

Since concrete cracks when subjected to relatively smal] 
tensile stress, a beam under negative moment essentially consists of 
the steel section and longitudinal reinforcement in the slab with the 
concrete acting as a transfer medium. For usual amounts of longitudinal 
Slab reinforcement the shear force transferred by the concrete slab 
through the shear connectors to the steel section is small in a negative 
moment region compared with that in a positive moment region. The be- 
havior of a composite beam under negative moment is essentially related 
to the characteristics of the steel section, the longitudinal slab re- 
inforcement and the shear connectors. In a negative moment region the 
lower flange is subjected to compressive stress, and local flange buck- 
ling may be a factor in the ultimate moment capacity. A number of com- 


(5) (6) 


posite beams tested under negative moment at the University of Alberta 
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and the University of Cambridge (10) (15) (21) 


failed in local flange 

buckling. It was observed in the University of AT pertaas tests that 
composite beams with a cover plate on the compression flange failed in 
lateral buckling because of increased resistance to local flange buck- 


ling. 


2.2 General Behavior of Composite Beams 
2.2.1 Positive Moment 


Since in a positive moment region the concrete slab is 
subjected to compressive stress, the composite section is similar to 
a steel section with a cover plate on the compression flange. Stress 
distributions for a composite section are illustrated in Figure 2.2 for 
three conditions, namely, prior to yielding, prior to strain-hardening 
and after strain-hardening. A typical moment-curvature relationship is 
Shown in Figure 2.3. 

In order to evaluate elastic stiffness for composite sections, 
it is customary to transform the concrete slab area into an equivalent 
area of steel by reducing the slab width by a factor, n = E/E: The 
bending moment derived from the idealized stress distribution shown in 
Figure 2.4 is defined as the theoretical ultimate moment capacity. At 
ultimate conditions the stress in the steel section is equal to the yield 
stress and the stress in the concrete slab is in terms of a Whitney 
stress block 22), The ultimate bending moment Mi for a beam with the 


neutral axis in the slab as in Figure 2.4(1) is therefore defined by 


Mi = A,oy(d/2 i t - a/2) Eid 
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where a = Ao, /0.85f (b. 


Mi for a beam with its neutral axis below the slab, as in 


Figure 2.4(b) is defined by 


My = Cie, + C.e, (Aig 
where C, = 0.85f A. 
and C,, = (Acoy - C )/2 


Effective interaction between concrete slab and steel sec- 
tion can be achieved by providing adequate shear connectors. Inadequate 
amounts or inadequate stiffness of shear connectors produces incomplete 
interaction. In the case of incomplete interaction two netural axes 
exist due to slippage as illustrated in Figure 2.5. The ultimate mo- 
ment capacity for a beam with inadequate shear connectors is obtained 
in a manner similar to that for a beam with adequate shear connectors. 
Since inadequate shear connectors cannot transmit the required force, 
the maximum compression force in the concrete is limited to the ulti- 
mate strength of the shear connectors, zQ.- Therefore a stress distri- 


bution at ultimate moment is as shown in Figure 2.6 and 
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The values of e, and e, in the above equations are dependent on the 
geometry of the beam section. 

Equations 2.2 and 2.3 form the basis for Figure 2.7 which 
represents the relationship between the ratio of the shear strength 


zQ/C; and moment capacity (23) (24) 


2.2.2 Negative Moment 


When a composite beam with longitudinal reinforcement and 
shear connectors in the slab is subjected to negative moment, such as 
in the vicinity of interior supports in a continuous span, the beam 
section may be considered to consist of the steel section and the lon- 
gitudinal reinforcement, since concrete cannot resist appreciable ten- 
Sile stress. The behavior is similar to that of a plain steel beam, if 
slippage does not occur between the slab and the steel section. Typi- 
cal stress distributions are illustrated in Figure 2.8 for conditions 
before yielding, before strain-hardening and after strain-hardening of 
the steel. Typical moment-curvature relationships obtained from Uni- 


(5) (6) 


versity of Alberta tests are shown in Figure 2.9. 
The simple plastic moment for negative moment conditions 


may be expressed as 


M =Ace, + Teo 2.4 
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where e and e, are defined in Figure 2.10. The simple plastic moment 
M, for a section, which has no longitudinal reinforcement, is the same 
as that for the steel section alone. 


(5) (6) 


Davison and Lever observed that the principal failure 


mode for composite beams subjected to negative bending was local flange 


8) 


buckling. Piepgrass observed lateral buckling failures in tests of 
composite beams with cover plates on the compression flanges. Figure 
2.11, based on University of Alberta tests, indicates that the ratio of 
the experimental ultimate moment to the theoretical ultimate moment 


decreases with an increase in the amount of longitudinal reinforcement. 


2.3 Properties of Material 
2.3.1 Steel 


Mechanical properties of steel significantly affect the 
behavior of composite beams. The typical stress-strain curve for a 
structural steel in tension consists of an elastic range, a plastic 
range and a strain-hardening range as shown in Figure 2.12(a) 629), 
The stress-strain relationship has been idealized by bilinear or multi- 
linear functions as illustrated in Figure 2.12(b). Function (1)»which 
idealizes the material as elastic and perfectly plastic, is employed 
in the simple plastic theory. Functions (2) and (3) express the stress- 
strain relationship by means of multilinear functions known as elasto- 
plastic strain-hardened functions. Function (2) is widely employed 


for structural analysis; however, since the tangent modulus after 


strain-hardening is defined as the initial tangent modulus at a 
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strain Est? Significant differences in strain occur between the approxi- 
mation and actual condition at higher stresses. In order to reduce the 
difference Lay and smith (26) proposed function (3) which approximates 
the strain hardening range by two straight lines. 

For purposes of the present investigation the stress-strain 
relationship after strain-hardening is idealized by function (4), which 
is a second order curve defined by three conditions; namely, the coor- 
dinates at the onset of strain-hardening (esp soy) the coordinates at 
ultimate stress (ce 20,) and the slope at the onset of strain-hardening 


Est: The stress-strain relationship may be expressed as follows: 


o = Ee for e<ey, 2.5(a) 
Oo = °y for eV act 2.5(b) 
o=atvb + Ce for ete 255 (c) 
Pons 2a = 
where i o 2E 40 (oF eu 
hE mere ay 3 (e - €.,)} 
u y Ste u st 
b = (oy - a)(o a cece) 


c= 2B, (oy - a) 


For example G40.12 structural steel with oy = 44ksi, ones 68ksi, Ent = 750ksi, 
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13 
Ec, = 0.0146 and oe 0.20, Equation 2.5(c) is defined by 


Oe= pol sot 7-40, 00tlS/20e 


The same relationship is assumed in compression. 
2.3.2 Concrete 


The stress-strain relationship for concrete, which has been 
the subject of extensive research related to the development of ulti- 
mate strength of concrete members, consists of an ascending portion 
and a descending portion as shown in Figure 2.13. 

In order to define the stress-strain relationship, various 
approximations have been proposed. noaneetadeabe proposed a parabolic 
equation for the ascending portion and a linear equation for the des- 


cending portion. 


me &c)\2 
f =f! (Per - Ee for ascending portion 2.6(a) 
Cc (e; Eo EG 
O.1be. 

f= F20/1.0.--———— for descending portion 2.6(b) 
e -€ 

cu Cc 

(28) 


Smith and Young proposed a single exponential function in the form 
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Where K is a constant and m is an exponent related to the concrete 


strength. A simple expression for the stress-strain relationship pro- 


posed by Desayi and Keishnans2) has the form 
Eve 
Ss Lt ae L260 
Cre at ele, 


where E. is the initial tangent modulus such that E. = ef feo: This 
expression is well fitted to the ascending portion as well as descending 
portion. 

Ultimate concrete strain based on results of axial compres- 
sion tests of plain concrete and beam tests of reinforced concrete has 
been discussed by many researchers. Ehenband reported an ultimate 
concrete strain of 0.0036 for cube strengths ranging from 3,000 to 7,000 


31) 


psi. Billet and nonleton, observed ultimate strains in beam tests 


varying from 0.0028 to 0.0040 with an average of 0.0038. Ros 34) | Brand- 


(34) 


EEO, Jensen and Hognestad et a1 (35) found that ultimate strain 


varied with concrete strength and proposed 


2 Rees 6 
eo 7 0-004 - £1/(6.5 x 108) 2.9 


from axial compression tests. A constant value Saree 0.003 is given 


in ACI Standard 318-71 36) (37) (38) for design purposes and is considered 
conservative. 
In contrast to the ultimate strain, the strain Eo at maximum 
(39) 


stress appears to be independent of concrete strength. Liebenberg 


observed values of Ec in the order of 0.002 in compression tests. 
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2.3.3 Shear Connectors 


Shear connectors are provided for two purposes in composite 
beams, i.e., transmission of horizontal shear between the concrete slab 
and the steel section and prevention of uplift of the concrete slab. 
Various types of shear connectors which fulfill these functions have 
been developed. Some types are rigid and give expectations of complete 
interaction. However, widely used types such as sarfh and channel con- 
nectors are flexible and permit slippage between the concrete slab and 
steel section. This slippage is due to deformation of the connectors 
and deformation of the concrete in the region of the connectors. 

Push-out tests have been widely used to determine the shear 
Strength and load-slip characteristics of shear connectors. The strength 
and load-slip relationship for connectors may be correlated with concrete 
strength. However, if load is expressed as a function of the static 
ultimate load, load-slip curves appear to be independent of concrete 


Strength as shown in Figure 2.14(40) 


(41) 


Yam and Chapman represented the load-slip curve by the 


exponential function 


q. = a {1 - exp(-ba,}} 2.10 


where q, is the ratio of shear force Q to ultimate strength QW? and 
a and b are constants. By choosing two points on the exponential curve 


so that Le = 2A., as in Figure 2.14, the constants may be defined as 
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In order to investigate the behavior of shear connectors 
in a negative moment region, push-out tests on studs set in reinforced 
slabs as shown in Figure 2.15 have been conducted by Johnson and Van 


(10) (13) (14) | Ultimate push-out strengths obtained varied between 


Dalen 
64 and 108 percent of the strengths obtained in conventional push-out 
tests, and slip at the failure load was approximately double that at 

80 percent of the failure load, although the load-slip curves were 
widely scattered. The behavior of studs was found to mainly depend 

on the location of cracks on the slab surface which influence the load- 
Slip relationship and the tensile force in the stud, resulting from the 


separation of the steel and concrete sections which influence the ulti- 


mate strength of shear connectors. 


2.4 Flange Local Buckling 
2.4.1 Previous Research 


Research on inelastic flange local buckling of plain steel 


beams has been conducted by various investigators, including Haaijer (42) (43), 


Lay (44) and Climenhaga $25) (21) | All these studies began with a consi- 


deration of buckling of an isolated plate in a region of yielding. 


(42) analyzed local buckling by plate theory. The 


Haaijer 
differential equation for an orthogonal anisotropic plate subjected to 


edge compression, as in Figure 2.16, has the form 
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where 2H Dyy Dix oS 4G,. The moduli Ds Ds Dy 


gested from the incremental theory of plasticity taking the second 


and G, were sug- 


invariant of the deviatoric stress tensor as the loading function. 
These moduli are functions of stress and strain, and Haaijer selected 
for his analysis the following values: D. = 3,200 ksi, DY = 32,800 ksi, 
Dyy = Dix = 8,100 ksi and G, = 2,400 ksi. 

tay es) analyzed local buckling by employing torsional 


theory as shown in Figure 2.17(a). The torsional buckling equation is 


in the form of 


ae) 


a46 d26 z doe 
Fey pie > BEY Geer MG ay pepe 2.12 


Lay evaluated G. on the basis of slip band theory which recognizes the 
basic discrimination of the behavior of steel in the region of yielding. 
This characteristic is not considered in the incremental theory. The 


ratio of inelastic to elastic shear moduli G,/G, was determined as 


ca 2 
hee airy) 68" a 


where a is the angle of slip plane as in Figure 1,17(b). Lay proposed 


the angle of slip plane is 45°, and Equation 2.13 becomes 
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The warping constant be and spring constant due to web buckling k were 


estimated as 


7 (b3t3 
Make ANS Pre | gue 
G.w3 
pues 2.16 
3(d - 2t) 


Using these values for Gs Je and k, the critical buckling load Eee 


is determined as 
mul nr) 2 L )2 


By differentiating the critical load, layer with respect to length L, 


the minimum buckling length is obtained as 


Since the local buckling shape is one full wave, the ratio L/b is ob- 


tained from Equation 2.18 as 
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There is some question as to the validity of this expression, since it 
indicates that the minimum buckling length increases as the web thick- 


ness decreases, which is contrary to what might be expected. 


(42) (44) 


Based on their respective analyses, Haaijer and Lay 
proposed maximum flange width-thickness ratios required to prevent local 
buckling. 


Haaijer proposed 


b__ $8.16 for 6 = 0.0 
2 ate for 8 = 0.1 ae 
where 6 is a constant for web restraint. 

Lay proposed 

Dawe 

aE 7 1.026VG,7o,, 2.21 


In his study of local buckling of composite beams, Climen- 


(15) (21) derived a relationship between the elastic-inelastic modu- 


haga 
lar ratio, h = E/E,» and the ratio of ultimate stress to yield stress, 


O,/o > on the basis of tests. He proposed 


ine) 


b= 0.348/ve, (3.18 - 07a) 
0.348/ Ey oh es Oy/oy Lee’. 


An approximate method for obtaining the moment-rotation 
relationship after local buckling was proposed by Climenhaga‘?>), based 
on yield line theory developed for limit analysis of concrete slabs. 


A local buckling wave was first assumed as a series of straight lines 
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derived from yield lines as shown in Figure 2.18, and the web deforma- 
tion was derived from a geometrical relationship for flange deformation. 
The moment-rotation relationship was then obtained by applying the princi- 
ple of virtual work. Subsequent composite beam test results showed 
reasonable agreement, although the theory is quite conservative for 


stocky sections. 
2.4.2 Proposed Analysis 


In the past considerable attention has been concentrated 
on establishing flange proportions for plain steel beams which guaran- 
tee that local flange buckling does not occur prior to the onset of 
Strain-hardening. It is necessary in the present study to develop a 
moment curvature relationship which extends into the strain-hardening 
range. 

The basic equation for torsional buckling of a plate sup- 


(45) 


ported by a spring at its centerline is expressed as 


dé de) _d {, jde = 4 | .2pde 
(El, I ; =(60§| ns ‘al ree eg 


in which the warping rigidity EI > St. Venant's torsional rigidity 
GJ, web restraint k and normal force P are all functions of z for a 
beam under varying moment. 

Figure 2.19 indicates the variation in tangent modulus EY 
with regard to strain. For a parabolic stress-strain relationship E. 
decreases with increase in strain as shown in Figure 2.19(b). The strains 


measured in the plastic region in a tension test represent over-all 
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strains, as would be measured over the gage length of a strain measuring 
device. Based on the theory of Luders lines,the plastic region is com- 
posed of an elastic zone, slip lines and a strain-hardened zone (42) (43) (46) (47) 
Therefore the tangent modulus in the plastic region is indeterminate. 
The value of the inelastic shear modulus G. in compression 


(42) (44) and others on 


has been proposed by Bijlaard’°). Haaijer » Lay 
the basis of deformation theory, incremental theory and slip band theory. 
Various expressions for G, are tabulated in Table 2.1. 

The finite difference equation equivalent of Equation 2.23 


at an arbitrary point i may be written as 
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The boundary conditions for a flange plate as shown in Figure 
2.20(a) may be considered as torsionally simple supports, that is, the 
rotation 6 and the second derivative of rotation 6" are equal to zero. 
However, the boundary conditions for a flange plate with multiple bear- 
ing stiffeners as in Figure 2.20(b) may be taken as a torsionally fixed 
Support, i.e., 6 =-6' = 0. 

The warping constant 7 for a thin rectangular section is 


)(150) and Saint-Venant's torsional constant J 


equal to b343/144 (149 
for a thin rectangular section is approximately bt°?/3.0 or more accu- 
rately bt3(1 - 0.63t/b)/3.06°2), The web restraint constant k accord- 


ing to Lay is G,w?/3.0(d 2 2t) ! 
2.5 Lateral Buckling 


Coos AeeLncroduction 


If the compression flange is reinforced with a cover plate, 


the resistance to local buckling is increased and lateral buckling may 
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occur prior to local buckling. prevorace 2) observed lateral buckling 
under negative moment in composite beams with a cover plate on the com- 
pression flange. The ultimate moment for the beams tested by Piepgrass 
was 3 to 7 percent higher than the plastic moment, and the curvature 

at lateral buckling was about five times the curvature at initial yield- 
ing. To simulate conditions in composite beams in negative bending, 


Climenhagatt>) 


tested steel beams with a cover plate on the tension flange 
in place of a concrete slab with longitudinal reinforcement. Lateral 
buckling occurred in beams where 5 x 3 RSJ steel sections were used. 

The lateral buckling configuration for a composite beam is 
illustrated in Figure 2.21. This configuration indicates that only the 
inverted tee section need be considered in a lateral buckling analysis. 
Elastic torsional buckling of a tee section with lateral restraint at 
the end of the stem has been analyzed for the case of axial load by 


Bleich (02) (53) on 


on the basis of thin-walled beam theory and Bulson 
the basis of plate theory. 

Lateral buckling of an inverted tee section in a composite 
beam due to axial force and bending may occur in the elastic-plastic 
range. Since plate theory is complicated in a plastic region, thin- 
walled beam theory is used in the elastic-plastic range. General equa- 


(49) (50) 


tions for thin-walled beams have been developed by Vlasov » Oden 


and others. Vlasov's theory, later modified by Rajasekaran and Varo) eet 
is applied herein. Based on the assumption that a cross section remains 
plane, a lateral configuration for an inverted tee section is idealized 


in Figure 2.21(b). 
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2.5.2 Lateral Buckling Analysis 


The stress distribution in a composite section under negative 
moment is shown in Figure 2.22(a). It is assumed that concrete can not 
transfer tensile stress, but the longitudinal slab reinforcement is effec- 
tive in tension. Isolating the inverted tee section, when lateral buck- 
ling occurs, an external force M9 must be applied along the end of the 
stem to balance the forces in the reinforcement and the upper flange 
as shown in Figure 2.22(b). The stress resultant in the inverted tee 
section may be separated into components due to normal force N, and due 
to bending moment M for the principal coordinates shown in Figure 2.22(c). 
When the section rotates about the hinge C, i.e., when the section 
buckles, additional indeterminate forces denoted as qx and qV in the 
x and y directions in Figure 2.22(b) occur at the hinge location C, be- 
cause of restraint in the x and y directions. 

For the above-mentioned conditions, buckling equations 
can be obtained from the general equations derived in Appendix A as 


follows: 


CE eel (Ee oO lias (Melati + (M)6")* = ay 2.25(a) 
(EI nt)" - {Nz(n' - a,8')}' + Tey, = ay 2.25(b) 


y 


(E1,9")" - (606')' = {Nz(~ayn' + ayé')}" + (ME")' - (M0")! 


are 


- a.) 2.25(c) 
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25 
where Cy and Cy are x and y coordinates of the hinge location, and ay 
and ay are coordinates of the shear center, which is located at the junc- 
tion of the web and the bottom flange. 
The coordinates Cy ays Xy and WE are zero, since the sec- 
tion is symmetrical about the y axis. Displacements - and Ne are 


zero iN the x and y directions at the location C. Therefore displace- 


ments € and n become functions of 6 as defined by 


ee ote 2.26(a) 


SNS 2.26(b) 


The coordinates Avs Cys Xp and w, are zero as shown in 
Figure 2.22(b). Substituting Equations 2.26 into Equations 2.25, the 


lateral buckling equation can be expressed as a function of 6 by 


(Elon) 2 - (GJe')' = Oe al) 


where ye = ie + (cy ~ ay 
to the hinge location C, and M, is defined as 


ae which is a warping constant with respect 


3 2 3 me 
M, Joss’ tely cy) }dA 2.28 


The local buckling equation, Equation 2.23, is developed as a par- 


ticular case of Equation 2.27 by setting C=ay and a constant stress 05» 
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and by adding a term for torsional restraint by web. The hinge loca- 
tion for local flange buckling is assumed to be at the junction of the 
lower flange and the web. 

The warping constant ae with respect to hinge location C, 
shown in Figure 2.22(b), is given by b3(d., + t/2)*/12 from thin-walled 
beam theory. The warping constant Lee for an inverted-tee section 


with a cover plate on the flange is evaluated as 


= 3 2 3 2 
Lgepy (4, + t/2)7t Hibs + tt 5/2) \/12 2.29 


The torsional constant J for a flange with a cover plate 


is defined as 


J = {1 - 0.63(t + t)/b, tb, (t + t,)°/3.0 + 


{1 - 0.63t, /(b, - b )} (by 


- b )t1/3.0 230 
where b,, and b, are the smaller and the larger of the widths of the 
flange and cover plate, and ty is the thickness of the plate corres- 
ponding to bv 

For the case where the difference between b, and bY is 
small, the second term may be neglected. For example, the ratio of 
the second term to first term is 2.7 percent for Piepgrass' test beams. 


The torsional constant may be expressed for such cases as 


J = {1 - 0.63(t + t)/b, b(t + t,)°}/3.0 
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FIGURE 2.20 BOUNDARY CONDITIONS FOR LOCAL BUCKLING AT AN INTERIOR 
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CHAPTER ITI 


DEFORMATION OF COMPOSITE BEAMS 


3.1 Introduction 


Beam deformations are produced by bending and shear. De- 
formation in composite beams is affected by slip along the interface 
between the concrete slab and steel section. 

Bending deflections are dependent on the moment-curvature 
relationship. This relationship is obtained by direct integration for 
a linear elastic stress-strain relationship. For non-linear stress- 
strain relationships a numerical integration procedure may be used. 
Several such procedures have been proposed, and most of them are based 
on dividing a cross section into a number of segments. Numerical inte- 
gration herein is based on the assumption of a linear stress variation 
in a segment, This assumption will provide adequate accuracy for a 
small number of segments. Since the stress varies in two directions 
for a section with residual stresses, the section is divided into seg- 
ments in two directions. 

Bending deformations may be calculated by means of such 
methods as the conjugate beam method, the finite difference method, 
the Newmark's method and the finite element method. The finite dif- 
ference method modified herein has several advantages for analyzing 


inelastic beam deformation. The method does not require the solution 
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of a large number of equations in each iteration since the number of 
unknowns equals the number of redundants. Although the conventional 


(57) 


finite difference method produces second order errors 0(\7 » where 
A iS Spacing, this error may be reduced appreciably by modification of 
the load term. 

Shear deformation can be evaluated independent of bending 
deformation, since it does not affect bending deformation, and in a 
continuous beam can be evaluated independently in each span. One of 
the simplest methods for evaluation of shear deformation in an inelastic 
beam is based on the principle of virtual work. This method is employed 
in the present study. 

The basic equation for slip deformation has been derived by 


Newmark (28) , vam(41) 


and others. The equation derived by Yam can be ap- 
plied to the deformation in an inelastic beam. However, the solution of 
the equation requires considerable computation time, especially for con- 
tinuous beams. The solution is simplified if we assume that slip strain, 


i.e., the derivative of slip with respect to length, is constant through- 


out the span. 


3.2 Moment-Curvature Relationships 


The stress resultants for a section, i.e., normal force 


N, and bending moment he are defined by 


N= fyodA 3.1(a) 
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52 
where the origin of the coordinates is the neutral axis as shown in 
Figure 3.1(a). 

Integration of Equations 3.1(a) and (b) can be performed 
for a linear stress-strain relationship without difficulty. However, 
a numerical method must be applied to more complex stress-strain rela- 
tionships such as those for concrete and for steel after yielding. 
Conventionally the integration of Equations 3.1 has been performed by 


a summation procedure 


ee 
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JpodA = Zo «dA. 32a) 


=x 
I 


= fnoydA = nosy dA. 3.2(b) 


where GV; and A. are the stress, coordinate of the centroid of dA; and 
area, respectively, as shown in Figure 3.1(b). 

Equations 3.2 do not provide sufficient accuracy for certain 
stress distributions. Accuracy is improved if us is defined as the 
coordinate of the center of the gravity for the stress block suggested 
by prnerde and shown in Figure 3.1(c). 

It is also possible to perform numerical integration by 
Newmark 's meted spe Assuming a linear stress distribution in a seg- 
ment as shown in Figure 3.2(a), the integration of Equations 3.1(a) and 
(b) results in 
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For equal areas of segments, the stress resultants N, and MY can be 


simplified as 


_ dA 
N, = are + ic, O54) 3.4(a) 
M = Gh 50g + Ao -etou.) Y 3.4(b) 
X 6 -1 jt 


It may be necessary to expand these expressions for normal 
force and bending moment to two dimensions for sections with residual 
stress. Denoting the equal intervals in the x and y directions by 
hy and ry and the stress at a mesh point by o44 as shown in Figure 


3.2(b), Equations 3.4(a) and (b) become 
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3.3 Effect of Residual Stresses on Bending 


Residual stresses occur in rolled sections as a result of 
plastic deformation during the cooling process after rolling. Usually 


compression occurs at the tips of the flanges and tension occurs at the 


(61) (62)(63)_ 


junction of flange and web in a wide flange shape Idealized 


residual stress patterns have been proposed by various investigators. 
The two examples shown in Figures 3.3 have been employed for analysis 


(63) (64) 


of torsional column buckling by Galambos and Lee et al “ 


Any residual stress pattern must satisfy equilibrium condi- 


tions, therefore the following equations must be eaciss eda 


fyopdA = 0 3.6(a) 


3.6(b) 
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In Equation 3.6(d) relating to twisting moment, the terms 
represent warping torque, Saint Venant's torque and a stress resultant 
which is defined by Equation A.25. If residual stresses are assumed 
constant along the member, the warping torque must be zero. Therefore 


Equation 3.6(d) can be rewritten as 


+M 9' =Q 3.6(e) 


which shows that twisting moment due to we balances Saint Venant's 
torque due to residual shear stress. If the twist is zero, that is, 


6' = 0, Saint Venant's torque must be zero. 


3.4 Bending Deformation 


The differential equation for bending deformation of a 


beam subjected to distributed load q, is expressed as 

(Ely")" = q, ony, 
Equation 3.7 can be divided into two parts 

Ely" = -M 3.8(a) 


Mu = -a, 3.8(b) 
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The finite difference equations for Equations 3.8 are ex- 


pressed in matrix form as 
FET) (BI {y} = -{M} 3.9(a) 


LAM teak 3.9(b) 


For a simply supported beam, the matrix [B] for equal seg- 


ments A can be established as 


The inversion of the matrix [B] is of the form, 


)* = 7b, 4] 3.11 


-1 
where the elements of the n x n matrix [B] , bij are 
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By multiplying by [B]~’, Equations 3.9 are expressed as 


{y} = -[B] *(—er)~+4m) 3.13(a) 


{Mi = -(B]7*{q,} 3.13(b) 


where [EI] "+ is a diagonal matrix with diagonal elements 1/EI.. 

Conventionally Equation 3.13(b) is substituted directly 
into Equation 3.13(a); however this procedure may yield second order 
errors 0(\2). As shown in Figure 3.4, constant moment is assumed in 
each segment, rather than varying moment as produced by loads. However 
the moment which is introduced into Equation 3.13(a) must be related to 
the actual bending moment diagram. By assuming a linear moment varia- 
tion between pivotal points i and (i+1), an equivalent moment at the 
point i may be given by = (M_, + 4M. + Migy): This may be expressed 
in matrix form by introducing co-efficient matrix [C] as 


4 1 
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Since this procedure is also applied to a distributed load, the equi- 


valent distributed load Veg and bending moment Mog are defined as 
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= [C){a,} 3.15(a) 
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[C] {M} 3.15(b) 


By using Equations 3.15 and substituting Equation 3.13(b) into 3.13(a), 


the deflection for a beam subjected to distributed loads is obtained as 


v= BED EB ei} 3.16 


Since a concentrated load P may be expressed as a distributed 
load equal to P/A, the bending moment {M} for a beam subjected to ex- 
ternal moments eb as concentrated loads {P/A} and distributed loads 


{q,} is expressed as 
(M} = (M} - (7 ({P/at + (elta,}) 3.17 


The corresponding deflection is 


fy} = BEEN i [-tm,,3+ [B]*({P/a} + [cI {a}) | 3.18 


Once the deflection {y} and moment {M} are obtained, the 
slope and shear can be evaluated by differentiating {fy} and {M} res- 
pectively. The shear Via at the pivotal point (i+) is expressed in 
finite difference form as 
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The shear V. at the location i is obtained by adding to the shear force 
Uh pee the additional shear AV. produced between i and (its), resulting 


in 


r 
Mia, 7 Mal/ a oareteas + 9541) ae 


where the second term is the shear at i for linearly distributed loads 


between i and (i+1). Similarly the slope can be evaluated as 


OM. OM. 
- y.)/r + e+ vs | AYA) 


Since the conventional finite difference equation 
does not include the co-efficient matrix [C], it implies that a unit 
matrix [E] is the co-efficient matrix, which is equivalent to a rec- 


tangular integration. 


3.5 Shear Deformation 


Shear deformation is often neglected, since it is not usually 
significant when compared with bending deformation. However, it is sig- 
nificant for short beams and for beams with small shear resistance. 

A composite beam has relatively large bending stiffness, 


but may not have large shearing stiffness, since the web of the steel 


08.€ (igGe * {PSS + AGH ‘ ai" a Sikes 


abso! batudtyteth yirsenth 7 t 35 npode ots at met bnosee ont anit : 
26 betsulsve od nso sok ont utrat ite: (DT) bas t negwied 


int ; 
‘ ; : phaay) 
[S.€ at, + aC uw) = 48 
. wee: iad ar Pe paptt Fale Ieee 


nottsups sonovattrb = st tnt? Isnotdnevncs eft gonte 
tiny 6 Seid estfgmt IF [3] xtrtom tnetoltts-oo arly sbufont ton 2e0b 
~95" 6 oF tnofevitups et dotdw .xtvdem Ingfottte-0o oft ef iE} xd. 7 
Hotisigotat vslugnss 


norssmotsd weede 2.6 


vl fso2u ton 2f Jt gonke “botoetpan nasto ef hobsoarvortab were ay : 
“whe et t .revoWnH not danttorsb, -pntbnad dttw bavsqmos nodw neg tingle 
-sonetetasn ~serle (Teme attw ensod ot brs amsed or2_ vot sneathin 
cezenttige pnibried apret Ulovitifen ze mad 93? 209n02 A es if 

fata ont to don sdt sonte .2zenttiye. Jon ¥ ne, 
s220nt = overt she 


1 7 i 7 , rie 
: 7 ss ian i : 
fa a a eh 4) ace et : ano 7 Aye 


section provides the major shear resistance. The ratio of shear de- 
formation to bending deformation also depends on the boundary conditions, 
and the type of loading. Bending deformation varies significantly for 
different boundary conditions. Shear deformation, however, is the 

same for all boundary conditions except for those at a free end or a 
yielding support. Shear deformation at z for the beam shown in Figure 
3.5 may be expressed as K.yp z{1-a) where Kon is the shear correc- 


tion factor defined by 
N 2 
_A ad) 


where Ny is the first moment of the section A, shown in Figure 3.6(b) 
with respect to the neutral axis. The shear correction factor is 1.2 
for a rectangular section, 2.94 for a W12x36 section, 6.63 for a com- 
posite beam consisting of a W12x36 and a 5"x60" concrete slab, and 3.1 
for a composite beam consisting of a W12x36 section and 3.6 in? of 
longitudinal reinforcement in a negative moment region. Generally Koh 
is close to unity for wide-flange sections if only the web is considered 
to resist shear. 

Shear deformation should not be neglected for an indeter- 
minate structure. For a beam with fixed supports the ratio of shear 
deformation to bending deformation is approximately 4 times that for a 
simply supported beam. 

Elastic shear deformation can be obtained by the virtual 


work method. Shear deformation is defined as 
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where the virtual shear stress tT equals to Vz/V, and V and V are the 
virtual shear and real shear. Therefore shear deformation Sch is re- 


presented by 


: | vy mda 3.24 
g 7A 


fntyda in Equation 3.24 can be derived from the normal stress au, due 
to bending at any location in the beam, since shear stress t is given 


by 
00 
Ate | 2 ap 3.25 
A 


The shear strain y is obtained from the stress-strain relationship 

which is similar in shape to the normal stress-strain relationship. The 
integration may be performed numerically by the same method as that ex- 
plained in Section 3.2. The computation of shear deformation is rela- 
tively simple, once (frydA)/V is calculated at intervals throughout the 


Span. 
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Analytical studies of composite beams with incomplete in- 
teraction were first undertaken by rege Sah who assumed linear 
stress-strain relationships for steel and concrete and a linear force- 
deformation relationship for shear connectors. The basic equation for 


the compression force C in the slab was expressed as 


S 
Cc" + K (FM - PIG} =NG 3.26 


where s is spacing, kK. is a constant related to the shear connectors, 


and 


pstoeret er Eta a7 


where Vs is the distance between the centroids of the slab and steel 
section. 

In an analysis involving non-linear stress-strain rela- 
tionships and a non-linear force-slip relationship, Yam and Chapman ‘42) 


derived two first order simultaneous equations with two independent 


variables C and A. (Figure 3.7) as 
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aoa =e TM} 3.29(b) 


where a and b are constants defined in Section 2.3.3. Yam solved Equations 
3.29 by means of the predictor-corrector method which is used for numeri- 
cal integration. For a simply supported beam a value of slip is assumed. 
However, for an indeterminate structure the bending moment must also be 
assumed. Since two assumed values are required for continuous beams, 
more iterations must be carried out to obtain correct values. 

Yam's approach requires considerable computational time. 
However, computational time may be significantly reduced by making an 
assumption that the slip strain, which is the differentiation of slip, 
is constant along the shear span. The slip distribution based on this 
assumption is linear as shown in Figure 3.8. The conditions in a shear 
span are illustrated in Figure 3.9 with all shear connectors in the 
shear span acting to resist the concrete compression force. The slip 
Strain is defined as cate A./L,; where A. is obtained from the force- 
Slip relationship and L. is the shear span. There are two neutral axes 
in a section due to slip, and the distance between them is defined by 
Ey/ >> where ¢ is curvature. The moment-curvature relationship can be 
obtained for this condition by employing numerical integration as dis- 
cussed in Section 3.2. Presumably the moment-curvature relationship 
depends onthe shear span and shear connector spacing. For the case 
of linear stress-strain relationships for steel and concrete and a 
linear force-slip relationship for shear connectors, the equivalent moment 


of inertia, accounting for slip, is derived in Appendix B. 
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Expressions for force-slip relationships in negative moment 
regions have not been proposed. Davison’?) measured strains in longi- 
tudinal slab reinforcement, and evaluated the ratio of measured strains 
to calculated strain, which he defined as an interaction factor. The 
interaction factor for the beams tested by Davison was approximately 


equal to 0.6. 


3.7 Deformations in Continuous Beams 


Stiffness in a negative moment region is usually smaller 
than that in a positive moment region, since in a negative moment 
region the concrete slab cracks when subject to relatively small ten- 
sion stresses. Therefore stiffness is dependent upon the bending 
moment, resulting in a material non-linear problem. 

Several methods have been proposed for the analysis of con- 
tinuous beams; however all methods require an iteration procedure for 
a non-linear problem. Finite difference and flexibility procedures 
are employed in the following. The flexibility matrix for finite dif- 


ference can be readily obtained for simply-supported beams as discussed 
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in Section 3.3. Continuous beams may be considered as a series of simply 


Supported beams to which the flexibility method is applied. 
ae 55 is the deformation at the support i of the primary 


structure, and Sas is the deformation produced at i due to unit value of 


redundant Xj> where oF is also defined as the flexibility factor, the 


j 
displacement A. at the support i becomes 
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or {a } = {6} + LFD{X} 3.30(b) 


Generally, moments or reactions are selected as the redun- 
dants in continuous beams. If moments are taken as redundants, the 
primary structure consists of n simple beams for an n-span continuous 
beam. One of the advantages of selecting moment redundants is that 
less computation is required in Equation 3.18. For example, a continu- 
ous beam with each span divided into m segments provides n{B] with 
m x m elements for a beam where moments are selected as redundants, 
instead of a [B] matrix with n(m + 1) x n(m + 1) elements for a beam 
where reactions are selected as redundants. 

The first step in the analysis procedure for continuous 
beams is to assume redundants {X}. The second step is to find {6} 
and [F] from Equations 3.18 and 3.26. Then redundants {X} can be solved 
from Equation 3.30(b). The third step is to find the moment in the 
real structure. If bending moments are close enough to those previously 
assumed, then the assumed redundants {X} are correct. If not, success- 
ive iteration must be carried out. The it assumption of {X} 1s pro- 


vided by 
{X}, = (Xb, + o({X}5_, > {X}y_2) Seal 


When the factor a is greater than unity, it is the so-called over-relaxation 
factor for elastic analysis. When it is less than unity, it is the so- 


called under-relaxation factor for inelastic analysis. 
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Several analytical methods have been proposed for an in- 
elastic region. One of the methods is based on using the stiffness 
EI in Equation 3.18 as defined by the secant modulus in the moment- 
curvature relationship. This method requires the evaluation of Equation 
3.18 in every iteration. As illustrated in Figure 3.9 , the total cur- 
vature for moment M can be divided into elastic and plastic curvature. 
Instead of changing the stiffness, an additional external moment, cor- 
responding to plastic curvature is applied to the inelastic region. 
This additional external moment is given by EI¢,, where EI is elastic 
Stiffness and by is plastic curvature. 

Since shear deformation does not affect redundants, shear 
deformation is determined after bending deformation. A computer pro- 


gram for deformation analysis is presented in Appendix D. 
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FIGURE 3.1 CO-ORDINATE SYSTEMS FOR NUMERICAL INTEGRATION PROCEDURES 
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(b) COORDINATE SYSTEM FOR TWO DIMENSIONS 
FIGURE 3.2 CO-ORDINATE SYSTEM FOR PROPOSED NUMERICAL INTEGRATION PROCEDURE 
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ASSUMED MOMENT DIAGRAM FOR 
CONVENTIONAL FINITE DIFFERENCE 
METHOD 


ACTUAL MOMENT DIAGRAM 


FIGURE 3.4 ASSUMED BENDING MOMENT DIAGRAM FOR CONVENTIONAL FINITE 
DIFFERENCE METHOD 
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FIGURE 3.5 SHEAR DEFORMATION NOTATION 
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FIGURE 3.7 SLIPPAGE MECHANISM (AFTER YAM AND cHapman' 41) ) 
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CHAPTER IV 
TESTS OF CONTINUOUS COMPOSITE BEAMS 


4.1 Introduction 


The object of the tests conducted as part of the present 
study was to provide additional information on the behavior of continu- 
ous composite beams. Although a number of tests have been conducted in 
the past by Barnardioeoe. Daniel et a1 (18) and park (42), it was felt 
necessary to conduct further tests with particular emphasis on the ef- 
fects of longitudinal slab reinforcement on ultimate strength, failure 
modes and general behavior. Three two-span beams were tested under 


concentrated loads applied at the midspans. 


4.2 Test Program 
4.2.1 Design of Specimens 


The size of beam specimens was determined to a large extent 
by limitations imposed by laboratory facilities. W12 x 31 and W12 x 27 
steel sections were selected, since these sizes have been used in pre- 
vious tests at the University of Mibents soe 4 CSA Standard 516-1969!) 
limits width-thickness ratio of compression flanges and webs for sections 
used in plastic design, and the ratios for W12 x 31, W12 x 27 and W10 x 21 
steel sections are approximately equal to these limitations. G40.12 


steel was specified. 
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18) 
4.2.2 Fabrication 


All steel sections were delivered to the laboratory complete 
with shear connectors and bearing stiffeners. Additional 3'-0" lengths 
of the W12 x 31, W12 x 27 and W10 x 21 sections were supplied for ma- 
terial property tests. The steel surface at locations for electrical 
resistance strain gages was ground smooth. 1" x 1" styrofoam blocks 
were placed at these locations and were subsequently removed after the 
concrete hardened. 

The slab reinforcing bars were also ground smooth at the 
strain gage locations, and then tied together to forma grid. Styro- 
foam blocks were attached to the bars at gage locations. Plastic chairs 
held the reinforcement in place. Plate 4.1 shows details of the fabri- 
cated steel specimens, reinforcing grids and form work in position prior 
to casting. 

Concrete was mixed in a 9 cubic foot Erich mixer in the 
University of Alberta Structural Engineering Laboratory. The mix propor- 
tions are shown in Table 4.4. One beam specimen together with 8 test 
cylinders was cast from 4 batches of concrete. The concrete was vibrated 
into place by a mechanical vibrator and was finished by means of a wooden 
screed and steel trowel. The beams were moist-cured for 5 days and 
removed from the formwork to cure under laboratory conditions for periods 
ranging from 69 to 72 days before testing. 

After a beam was positioned in the testing apparatus and 
strain gages were connected to the indicators, the steel section was 


white-washed so that the development of yielding could be observed. 


ne ee, io 
; _ " ; a ; : 
mat ST, o: oan) Ce , m 
{ ¥ ia " i i i gs i] c : = ; 
| Ave to eee had | viv outset visi $.S.8 
. 1 t 


atelqmos yrotsiods! eft ot beveviteb ovo anottose fast2 TTA 
edjpnal “O-'€ IsnokstbbA .2tsnetths2 pnivssd bas etotsenn09 ssade dat 


~6m “oT bei fqque svaw enoftose IS x OW brs YS x SIW «f& x SI ott 40 
[sottto9Te or enotssso! 35 aostwe foat2 aft .2te0d ytveqotg isiie3 


aioold maotoryte “#I x “AI .Adoome bavotp e6w 29p6R ates conetst oe? - 
tteo0! seodt 36 baoslq ovew 
* benebysd st979009 


gnats p 


afd vette bevomer yIingupgedue ovew brs 20 


edt 2s dtoome bruoww o@fs o19w 2766 ontovotntea dele ont 
~orgsd | .bive 6 mot oF tedspod bei? nate brs .2notssoot spse ntsise 
axvtaits otteal? .enoftsool eps ds cved sid 0s boriastys S19" zd0fd msot 
-tdst od Yo 2tteteb awone I.6 otei4 .o26!q at SnomoertoTnrsy ot bisd 
sort nottteoq nt drow mo? brs abite prfovotnioy _2nembosqe Teste bats. 
spatt2s9 oF 

oft nt vextm dof1a toot ofdu> © 6 nt boxtm enw stavoned 3 
-rogorg xtm oaT .iotsvede! ontysentend [ swiourse sdvedlA to ystetevind 
test 8 itt “aridapos nomtosqe mead ond = .h.A gfdeT at wore 616 2nord 
botevdtv 2ew sisvanos afl .steron0o to esrlodad & mor? Jes 26 erabnt Iya 
naboow 5 to ensem yd barietnit 26w brs votsydiv feotnsoom 5 vd s06fq odnt 
bag eysb 2 vor -betusstetomovew ensod sft .teword’ Teste brs beerde 
zbot1aq rt enottibnos yrossveds! vebau sw os Aowniot ony mov? bevomet 
-pnttest aroted 2ysb SV oF ed mov? pnftpast 

bas avtsisqgs entteet old nt benofdteeq 26w mood 6 W9StA 

z6w notice2 fset2 anv ,erossotbnt ats oF besasano> svew 2eRsR nterwe 
veenntty coh nan aaeermne ?* snamgotaveb orft Jedd oz beresw-adirw 


rh) ren ole ; ess ’ onh-48 17) rheee 


78 


A 4" slab thickness and a 4'0" slab width were selected, 
Since these dimensions had been used in previous tests at the University 
of Alberta. A concrete strength of 4500 psi, which is a common strength 
used in actual structures, was selected. 

The amount of longitudinal slab reinforcement at the interior 
support significantly affects ultimate loads and failure modes. Local 
buckling may be the dominant failure mode in a beam with a relatively 
large amount of longitudinal reinforcement. On the other hand, crushing 
of concrete may be the dominant failure mode for a section with a smal] 
amount of longitudinal steel. The amount of longitudinal slab reinforce- 
ment as shown in Table 4.1 were selected to produce different failure 
conditions. The longitudinal bars were placed at mid-depth in the nega- 
tive moment regions. 

Shear connectors as shown in Figure 4.1 were provided ac- 
cording to provisions of CSA Standard S$16-1969, which are based on 
ultimate connector strength. The shear connectors were 3/4 inch diameter, 
344 inches long,headed stud connectors. 

Transverse reinforcement was introduced in order to control 
longitudinal cracking. samaae CU! encountered longitudinal cracking 
in beam specimens with transverse steel in the amount of 0.2 percent of 
the concrete area,which corresponds to temperature reinforcement require- 
ments. The present beams were reinforced transversely with #3 bars 
placed at 43," intervals throughout the length of the slab. This rein- 
forcement is equal to 0.67 percent of the concrete slab area or three times 
the nominal temperature reinforcement. 


Complete details of the test specimens are shown in Figure 4.1. 
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4.2.3 Test Equipment 


Test specimens were supported at mid-length on a hinge reac- 
tion unit anchored to a concrete pedestal and at each end by a roller 
unit fitted with a rocker plate assembly seated on a concrete pedestal , 
as shown in Figure 4.2. 

Two 100-ton Amsler hydraulic jacks were used as loading units 
centered on the longitudinal centerline of the specimen. Loads were 
transmitted to the concrete slab through 12" x 8" x 1's" spreader plates 
set in plaster of Paris to ensure uniform contact. 

Details of the lateral support system are shown in Figure 
4.3(a) and Plate 4.2(a). At each support, steel channels, attached to 
the web stiffeners and supported by vertical steel columns, resisted 
lateral rotation of the specimen as shown in Figure 4.3(b) and Plate 
4.2(b). The slab was supported laterally by rollers on vertical guides 


near midspans. 
4.2.4 Material Properties 


Results of tension tests performed on samples of #3, #4 
and #5 bars used as longitudinal and transverse reinforcement are shown 
in Table 4.2. Test coupons, two from each flange and three from the web, 
were cut from 3'-0" lengths of W12 x 31, W12 x 27 and W10 x 21 provided 
for this purpose. Results of coupon tests are shown in Table 4.3. 

Eight concrete cylinders were cast for each beam and cured 
under laboratory conditions. Five cylinders were tested in compression 
and three were subjected to a splitting test at the time of the test on 


the corresponding beam. Test results are shown in Table 4.4. 
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4.2.5 Instrumentation 


Measurements of strains in the steel sections, reinforcing 
bars and concrete slabs were obtained by means of KFC-5-C1-11 and SR4- 

A7 electrical resistance strain gages. Deflections, rotations and slip 
deformations were measured by mechanical means. 

Deflections were measured by means of a precise level fo- 
cused on graduated scales suspended from the bottom flange of the speci- 
men. Deflection measurements were taken at the centers of the two spans 
and at locations 30 inches from the interior support in the vicinity of 
the points of inflection. 

Rotations were measured by means of mechanical rotation meters 
at the beam ends and by means of sets of two dial gages at locations 
30 inches from the interior support. These mechanical rotation meters 
consisted of a level tube, an extension dial and two arms connected by 
a hinge as shown in Figure 4.4. A vertical steel bar was welded to the 
steel section and two dial gages 24 inches apart were mounted by magnetic 
bases on the interior support as shown in Figure 4.4. Rotations were 


determined from the relationship 


where Ad and Ag, are displacements measured by the dial gages and Le is 
the distance between the gages. 
Curvatures were measured at the interior support and the 


load points by means of electric resistance strain gages. For measuring 
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curvature, three gages were attached to the steel section, six gages were 
attached to the concrete slab in positive moment regions and three gages 
were attached to reinforcing bars in the negative moment region. 

Points of inflection were determined from measurements ob- 
tained from sets of two electric resistance strain gages. Since in- 
flection points were expected approximately 30 inches from the interior 
Support, gages were positioned at 24" and 36" from the support. Assuming 
that strains in the flange vary linearly, the inflection point can be 
simply located from the measurements obtained from these two gages. 

Once inflection points are located, reactions and bending moments are 
easily evaluated from an analysis of a statically determinate structure. 

Slip and slip strain measurements were obtained by means 
of dial gages and electric resistance strain gages, respectively. Six 
locations for these measurements were selected in the positive moment 
region. Slip strains are defined by the difference of strains in the 
steel section and concrete slab at the interface. Therefore two strain 
gages were required at each location. Slip was measured by means of 
dial gages attached to the steel sections by magnetic mounts, and react- 
ing against a steel angle attached to the concrete slab. 

Strains in the transverse reinforcing bars were measured 
at four locations. 

Complete instrumentation details are shown in Figure 4.4 and 


Plate 4.3. 


4.2.6 Testing Procedure 
Prior to testing, the electric resistance strain gages were 


connected through switch boxes to strain indicators. Immediately prior 
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to loading all gages were balanced for initial reading and mechanical 
gages were positioned securely. A load of 5 Kips was applied through 
each jack and then removed before initial readings were taken. Existing 
cracks in the slab were marked at this time. 

Load increments were different in each test. For steel 
strains less than the yield strain 10 Kipsor 20 Kips increments were 
applied. After yielding began, load increments were gradually reduced 
to provide sufficient values to plot an accurate load-deflection relation- 
ship. For each load, strain and deflection measurements were recorded 
and cracks were marked. At each load sufficient time was allowed for 
stabilization of deflection and rotation. After maximum loading the test 
was continued in order to obtain the falling portion of the load deflec- 
tion curve. When deflection approached the 5 inch maximum extension 
of the jacks, load was removed and additional 1 inch plates were posi- 
tioned under the jacks and load was re-applied. Beams were tested 
well beyond ultimate load conditions to the point where concrete crushing 
and/or significant flange buckling at the interior support occurred. 
Failure modes and crack patterns were photographed before the specimen 


was removed. 


4.3 Test Results 


4.3.1 Introduction 


All original data obtained from the beam tests is filed in 
the Department of Civil Engineering at the University of Alberta. Data 


presented herein is in tabular, graphical and photographic form. 
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4.3.2 General Behavior 


4.3.2.1 Beam CB1 

No visible cracks could be found in the concrete slab at the 
Start of the test. At a load of 30 Kips, cracks appeared on the top 
surface of the slab in the negative moment region. These cracks ex- 
tended to the bottom surface of the slab at a load of 50 Kips. Cracks 
appeared on the bottom surface in the positive moment regions at a load 
of 50 Kips. Spalling of the white-wash on the steel section began at a 
load of 80 Kips in the negative moment region and at a load of 100 Kips 
in the positive moment regions. 

At a load of 100 Kips it required approximately 5 minutes to 
establish equilibrium, and more than 15 minutes were required to sta- 
bilize the load at 120 Kips. Cracking progressively spread throughout 
the slab and crushing began at the top surface at a load of approxi- 
mately 131 Kips. At a load of 135 Kips it was still possible to stabi- 
1ize conditions in the beam, though crushing extended over 15 inches. 
The load was then increased to 137 Kips. After 15 minutes the slab was 
crushed completely, and the load could not be held. After crushing of 
the concrete, local flange buckles formed and the load decreased. 


Crushing and local flange buckling failures are shown in Plate 4.4. 
4.3.2.2 Beam CB2 


This beam was the first beam tested in this program. It 
was loaded initially without a bracing system. However, at a load of 
70 Kips the beam rotated transversely to such an extent that the test 
was terminated. After the bracing system was introduced, the beam 


was load to failure without any further difficulty. 


at de > : va 7 a 

. ) ii, 
| | var ; rae wa 
ue i o> 

te oe z e ae 

OW oe Eb ; 


at © me oa af 
sit t6 dele aseranos aft nt bnuot 2 biuos becccagitie a 
masts 14 ae 


i) 
& 
W 


6 t6 mispad nottose fate ont no desw-od-biw 9nd +6 ont tage 24h ae : 7 
autN.OOf to béol 6 3h bns noisy taamom ovtyapan ons tt ant 0&8 to bsof 
Be” a) | 
Lenotne, neon evttt20q any nt 
PT VaGre 


ot eatuntm ¢ vlstemixotaqs bovtupey, Sf eqrA OOF to bsof 6 JA 
' (ée 
-st2 03 battupay ovsw eotuntm ef ney svom bis .muradtltupe det idss2s 


a- 7 aye 9 
tuonpuorilt beatg2 yfevi2zestporg prftosid .2qht USE ts bsof ont ike 
13 


-ixonqgs to bsol 6 de Sostwe qos ond Js msped pf Fitaut brs dole ong 
3 

date oF afdizeog PMtte 25w It eqkX céf To beol 6 JA .2qta rel wags 

.eariont eI ysvo baba oniizuy> dovodd .mesd odd nit anor Hbn0> ext . 


saw defe att 2otunim ef westA .2eqra NEL of bezsorvont nent 26w bsof oat 
Pek Ue oy 
to prtdauya vettA «bled sd ton bios bsof snd bas Xfatotqmoa bereu12 ep 
ey eae ’ 


_baeserpeb bsol edd bas bemrot 2eftoud spnstt fsoof ,sd9T0N09 ods 
fp. 8 otsf4 at wore ois 2awwltst pnitdaud spnef? fsa0f bas ent dawns 


a 
” ' 


$63 msod $§.S.E.8 
a Y & 


’ 
+i .meypova 2etdt af bestest mood t2vtt ord 2sw meed ete . 


to bsof 6 ts _19vBWOH ‘nas 242 pniosid 6 Suortiw xitsra tar ee 
teed att dent Ins3xe 16 dave of Visevevensnd, badstor mad oft eqt 1 


wh tole 
msod sft Saati 7) deve a ante ® 8 aT bstentim 9 a " 


84 


The first cracks were observed at a load of 20 Kips in the 
negative moment region and at a load of 50 Kips in the positive moment 
regions. Local flange buckling initiated in the negative moment region 
at a load of 120 Kips and extended gradually with increasing load. Local 
buckling was completely formed at a load of 130 Kips. At this load 
crushing of concrete started at the load locations. Deflections were 
Stabilized at a load of 133 Kips, but could not be stabilized at 136 
Kips. The load decreased to 130 Kips at which point crushing was ini- 
tiated. During unloading, web buckling occurred between the load point 
and the interior support at a load of 126 Kips. Finally a section of the 
concrete slab at a load point spalled off at a load of 116 Kips. The 


failure modes are shown in Plate 4.5. 
4.3.2.3 Beam CB3 


A visible crack existed on the top slab surface near the 
interior support prior to loading. Tension cracks began to appear in the 
negative moment region at a load of 10 Kips, and in the positive moment 
regions at 50 Kips. Local buckling initiated at a load of 87 Kips and 
Slowly developed as the load increased. At a load of 92 Kips crushing 
of concrete occurred at a load point and at a load of 94 Kips the con- 
crete was completely crushed. The flange buckles were completely formed 
at the interior support as the deflection increased. Failure modes are 


shown in Plate 4.6. 
4.3.3  Load-Deflection and Load-Rotation Relationships 


Load-deflection relationships are shown in Figures 4.5. 


Deflections were measured at midspans and at locations 30 inches from 
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85 
the interior support. Rotations at the exterior support and at loca- 
tions 30 inches from the interior support are plotted against load in 


Figures4.6. 


4.3.4 Load-Reaction and Load-Moment Relationships 


Inflection points were determined from strain data obtained 
at locations 24 inches and 36 inches from the interior support. Reac- 
tions were then calculated on the basis of a statically determinate 
structural analysis. Resulting load-reaction relationships are shown 
in Figures 4.7. Bending moments at the load locations and at the interior 
Support were obtained from reactions and loads. Load-moment relation- 


Ships are shown in Figures 4.8. 
4.3.5 Moment-Curvature Relationships 


Moment-curvature relationships at the midspans and at the 
interior support are shown in Figures 4.9. The maximum strain values 
obtained from electric resistance strain gages were 2 to 3 percent. This 
magnitude of strain was sufficient to establish moment-curvature relation- 


ships. 
4.3.6 Slip Deformation and Slip Strains 


Load-slip relationships are shown in Figures 4.10. Distri- 


bution of slip and slip strain are shown in Figures 4.11 and 4.12. 


4.3.7. Strain Distribution Across Slab Width 


Strains were measured at the slab edges andthe longitudinal 
slab centerline. Figure 4.13 shows the ratio of slab edge strain to slab 


centerline strain at the load points for various load values. 
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4.3.8 Transverse Strains 


Strains were measured in the transverse reinforcing bars at 
the longitudinal centerline and at various locations in the span. Load- 


strain relationships were shown in Figures 4.14. 


433.9 Crack Patterns 


Crack patterns in the concrete slabs are shown in the Figures 


4.4 Discussion of Test Results 


4.4.1 General Behavior 


Transverse rotation may occur as a result of warping of the 
slab or steel section prior to testing. Transverse rotation due to load 
was effectively resisted by the bracing system and by the reaction system 
described in Section 4.2.3. 

Initial cracks were observed only in Beam CB3. This crack 
may have been due to shrinkage, or, since it was located over the interior 
Support, it may have occurred during positioning in the test apparatus. 

The first transverse cracks appeared in the negative moment region. Later 
transverse cracks appeared on the bottom surface in the positive moment re- 
gions. Diagonal cracks developed at the approach of failure. Longitudinal 
splitting was effectively prevented by means of the transverse reinforcement. 
Although a few longitudinal cracks appeared near the load locations before 
failure, they did not appear to affect the ultimate moment capacity of the 
beam. 

Failure of Beam CB1 was a concrete crushing failure. Local 


buckling followed as a result of reduced stiffness in the positive moment 
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87 
regions. The amount of longitudinal reinforcement in Beam CB1 was 1.6 in? 
or 17 percent of the area of the steel section, which was the smallest amount 
in the three beams. Failure of Beam CB2 occurred as a result of local buck- 
ling followed by concrete crushing. The area of longitudinal reinforcement 
was 3.1 in*, i.e., 38.9 percent of the area of the steel section which was 
the highest amount in the three beams. Failure of Beam CB3 resulted from 
concrete crushing followed by local buckling. Although local buckling ini- 
tiated prior to crushing of concrete, it did not produce a sudden failure. 
The area of longitudinal reinforcement in Beam CB3 was 32.3 percent of the 


area of the steel section. 
4.4.2 Load-Deflection and Load-Rotation Relationships 


Load-deflection and load-rotation relationships have similar 
characteristics as indicated in Figures 4.5 and 4.6. They consist of three 
ranges, i.e., elastic, plastic and unloading ranges. The deflections at the 
maximum load were 2.6, 1.8 and 2.1 inches for Beams CB1, CB2 and CB3, respec- 
tively. These values indicate an increase in deflection with decrease in the 
ratio of the area of longitudinal reinforcement to the area of the steel sec- 
tion. 

Rotation measurements at the ultimate load were not dependable. 
Since the rotation bars were welded to the web of the steel section, deforma- 
tion of the web affected the rotation measurement. End rotations at ultimate 
load were approximately 5 times the rotations at initial yielding. The un- 
loading portions of the load-rotation curves are different to those obtained 
for isolated simple beams under negative bending due to the fact that con- 


tinuous beams have greater ductility than isolated simple beams. 
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4.4.3 Load-Reaction and Load-Moment Relationships 


Reactions were determined from the location of points of 
inflection. The moment diagram is linear for concentrated loads. If 
the concrete slab does not have tensile cracks, the strain can be assumed 
to vary linearly along the span and the location of a point of inflection 
may be determined by linear interpolation between strains measured on 
both sides of the point of inflection. If the concrete slab has tensile 
cracks, the location of a point of inflection is related to the stiffness 
and the distance from the neutral axis to the gage location in positive 
and negative moment regions. However, the effect of cracking to moments 
is not significant. Strains 24 and 36 inches from the interior support 
did not exceed yield strain values in any of the beams. The points of 
inflection in each span were almost equidistant from the interior sup- 
port and therefore the two exterior reactions and the moments at the 
load positions were almost equal in the elastic range. However, in the 
plastic and unloading portions they differed slightly as shown in Figures 
4.7 and 4.8 probably due to imperfections in the beams affecting symme- 
try. 

Bending moments are proportional to the applied load in the 
elastic range. Figures 4.8 show slight moment redistribution in Beams 
CB1 and CB2 for loads greater than 90 Kips and significant redistribution 


in Beam CB3 for loads greater than 50 Kips. 
4.4.4 Moment-Curvature Relationships 


Curvature values were obtained from strains at the top and 
bottom of the steel sections. The range of curvatures which could be 


measured by this means was limited because the strain gages were not 
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effective beyond strains of 2 to 3 percent. Therefore curvatures were 
not obtained for the unloading range. 

Although moment-curvature relationships for Beam CB1 indicate 
that stiffness in the positive moment regions was greater than in the 
negative moment regions for elastic conditions, inelastic stiffness was 
approximately the same in both regions. Therefore the ratio of inelastic 
to elastic stiffness in the negative moment region is greater than that 
in the positive moment regions. This may be related to the descending 
portion of the stress-strain relationship for concrete. This effect is 


apparent in all beams tested. 
4.4.5 Slip and Slip-Strain 


Maximum slip occurred at locations near the load, as shown 
in Figures 4.10. This behavior is similar to that observed in simple 


(20) (40) Since shear forces are in opposite directions on 


span beams 
opposite sides of the load, slip deformations reverse at the load loca- 
tion. Load-slip relationships presented in Figures 4.9 include the un- 
loading portion, which is not defined in the pushout test. 


As shown in Figures 4.11, there is no particular trend in 


the slip-strain distribution. 
4.4.6 Strains in Transverse Reinforcement 


Strains were measured in the transverse reinforcing bars 
which were located 3 inches from the bottom of the slab in positive mo- 
ment regions and 2% inches from the bottom of the slab in negative mo- 
ment regions. Figures 4.14 indicate that strains in the transverse re- 
inforcement exceeded yield strain in positive moment regions but not in 


the negative moment region. 
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4.4.7 Effective Slab Width 


The effective slab width may be evaluated from the stress 
distribution. The ratio of strain at the edge of the slab to that at the 
center, y.> was approximately 70 percent as shown in Figure 4.13. The 
theoretical strain distribution over the slab width is an exponential 


(51). 


function Assuming the function as a parabola, the effective slab 


width, ba» may be expressed as 
Mahe s ey a 4.1 
e CaaS Yeit'c 


For the beams tested, b = 6.5 in., DA = 48 in. and eae Osi 
Therefore the effective slab width is 39.7 inches. This results in 
(b, - b)/(b. - b) equal to 0.81 or in (b. - b/&equal to 0.23. 

According to CSA Standard s16(7) the effective projection 
of the slab beyond the flange of the steel section is considered to be the 
smaller of: (i) one-sixth of the beam span and (ii) twelve times the 
Slab thickness for a slab which is not supported along its edges, and 
the smaller of: (i) one-fourth of the beam span and (ii) sixteen times 
the slab thickness for a slab which is supported along its edges. The 
effective design width, therefore, is 30.5 inches for the former case and 
40.5 inches for the latter case. The effective slab width for the test 
beams with unsupported slab edges was therefore greater than the Code 


value and was close to the Code value for a beam with supported slab edges. 
4.4.8 Crack Patterns 


Three crack patterns were developed in the slab, ij.e., 
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transverse cracking, longitudinal cracking and herringbone cracking 

as shown in Figures 4.15. Transverse cracks were developed in the nega- 
tive moment regions due to tensile stress at relatively small loads. 
Longitudinal cracks at the load points were produced along the centerline 
of the beam at relatively large loads. Herringbone cracks were developed 
at large loads and the direction of the cracks was approximately 45° to 
the beam centerline near the ends of the beams. The angle decreased as 
the load point was approached. These two crack patterns, i.e., longi- 
tudinal and herringbone, are caused by the stresses produced around shear 


(65,66) 


connectors and by the relatively large tensile strain as a result 


of Poisson's ratio effect. 
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TABLE 4.2 PROPERTIES OF SLAB REINFORCEMENT 


BAR SIZE tae pte 


ULTIMATE STRESS 
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TABLE 4.3. PROPERTIES OF STRUCTURAL 


STEEL 


STRAIN-HARDENING 
MODULUS 
ksi 


1330 
1250 
1330 
1360 
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1140 
1470 
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1160 
1000 
1130 
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1000 
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850 
700 
750 
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TABLE 4.4 


AGE AT 
BEAM TEST 
days 


if 
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CB3 72 


* MIX PROPORTIONS 
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WATER 
SAND 


COARSE AGGREGATE 


PROPERTIES OF CONCRETE* 


COMPRESSIVE 
SHERESS 
psi 
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6366. 
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a33R1DSe 
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344 lbs. 
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(b) REACTION SYSTEM 
FIGURE 4.3 SUPPORT DETAILS 
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PLATE 4.1 | FORMWORK AND REINFORCEMENT DETAILS 
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(b) SUPPORT SYSTEM 


PLATE 4.2 LATERAL BRACING AND SUPPORT SYSTEM 
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PLATE 4.3 INSTRUMENTATION 
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PLATE 4.4 FAILURE OF BEAM CB1 


bel 


: 
i 
; “Eee 
; Ore 
erie 


ie ay” 
aiaulaatl 
rT mt hoe y, - 

i ; , 


- 
: 


AS 4 etl] ¥: a | 
ae. 

‘ash 

4 


4 


ae 
: 7 te 6) 
ms 


£40 MARE 40 JAUIIAR =$.8 STAN 


135 


BEAM CB-2. 
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PLATE 4.5 FAILURE OF BEAM CB2 
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PLATE 4.6 FAILURE OF BEAM CB3 
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CHAPTER V 
BEHAVIOR OF CONTINUOUS BEAMS 


5.1 Introduction 


The present chapter compares results of the proposed analysis 
with experimental data from tests conducted at the University of Alberta, 
the University of Cambridge and Lehigh University. 

Deformations in continuous composite beams have received 
relatively little attention in the past. Deformations have been mainly 


(8) 


analysed assuming elastic properties. Barnard discussed deformation 
and moment-curvature relationships without consideration of moment redis- 
tribution. In the present chapter deformation behavior is discussed for 
both elastic and plastic ranges. The analysis includes the effect of 
Slip and shear as well as bending. Moment redistribution, which occurs 
in the process of the formation of a mechanism, is taken into account. 
Crushing of concrete in the positive moment region and local 
flange buckling in the negative moment region are the major failure modes. 


Shear connector failure and longitudinal splitting can be prevented by 


proper design provisions. 


5.2 Local Flange Buckling 


5.2.1 Presentation of Data 


Havisohis and Riea beams were composed of W12 x 27, 
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W12 x 31 and W12 x 36 sections and concrete slabs with various amounts 
of longitudinal slab reinforcement. All beams failed in local flange 
buckling after ultimate moment capacity was achieved. Dimensions of test 
beams are shown in Figure 5.1. 

Theoretical values shown in Table 5.1 are based on Lay's 
and Bijlaard's proposed shear modulus values. Figures 5.2(a), (b) and 
(c) compare test and analytical ultimate moment values with plastic moment 
values. They also illustrate the effect of longitudinal slab reinforce- 
ment on theratio of test ultimate moment to the plastic moment values. 
Table 5.2 compares test and theoretical values of curvature at the onset 
of local flange buckling. Theoretical values are based on Lay's and Bij- 
laard proposed shear modulus. The relationships between curvature 
based on Lay's proposed shear modulus and the amount of longitudinal rein- 
forcement are shown in Figure 5.3. 

To simulate the effect of longitudinal slab reinforcement in 
a composite beam under negative moment, Clinenhaga (15) tested beams 
consisting of a steel section and a coverplate welded on the tension flange. 
Some beams contained longitudinal stiffeners welded on the web as indi- 
cated in Figure 5.4. Test and analytical results are tabulated in Table 


Dade 
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For beams tested at the University of Alberta, experimental 
and theoretical values for local buckling moments are very similar. Lay's 
inelastic shear modulus results in higher moment values for stockier 


sections such as W12 x 31 and W12 x 36. Generally theoretical values 
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based on Bijlaard's shear modulus are higher than those based on Lay's shear 
modulus in the vicinity of the onset of the strain-hardening, and tend to 
decrease with increasing strain. Composite beams with W12 x 27 sections 
appear to buckle shortly after the onset of strain-hardening. 

The difference between experimental moment values and theore- 
tical moment values based on Lay's shear modulus for beams tested at the 
University of Alberta is within 15 per cent. The difference is large 
for beams with W12 x 27 and W12 x 31 sections and lower amounts of rein- 
forcement, and also for beams with a W12 x 36 section with larger amounts 
of reinforcement as shown in Figures 5.2. The figures also indicate that 
theoretical values increase in proportion with the simple plastic moment 
values for increasing amounts of reinforcement. However experimental 
values do not increase as much with increased amounts of reinforcement. 
Values of MM, indicate this behavior as shown in the same figures. 
Although for beams with W12 x 36 sections experimental and theoretical 
values are fairly close, they still exhibit the same above-mentioned 
trend. This may be due to underestimating the web restraint for small 
amounts of longitudinal reinforcement, or due to the effect of lateral 
displacement for large amounts of longitudinal reinforcement. However, 
since the difference is not significantly large, the proposed analysis 
is considered acceptable for predicting local flange buckling. 

Figures 5.2(a), (b) and (c) show that the ratio of ultimate 
moment to the simple plastic moment, M,/M,» for beams with W12 x 36 sec- 
tions is larger than that ratio for beams with Wl2 x 31 and W12 x 27 sec- 
tions. This is presumably due to the large flange torsional resistance 


of W12 x 36 sections. The width-thickness ratio of a W12 x 36 section 
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is 6.1 which is considerably larger than that specified in CSA Standard 
$16-1969(1), 

Table 5.3 shows good agreement between theoretical and experi- 
mental ultimate moment values for Climenhaga's beams, except for beams 
with 5x3 RSJ sections and for beams with longitudinal stiffeners on the 
webs. The average ratio of test to theoretical moment values is 1.11. 

The shallowest sections (5 x 3 RSJ) exhibit the largest ratio of test 

to theoretical values. The presence of longitudinal web stiffeners increases 
the BRST ad) approximately 10 percent for 16 x 5%; UB26 sections 

and by approximately 20 percent for an 8 x 5% UB17 section. A comparison 

of experimental and theoretical moment values for 8'-4" and 12'-4" spans 
indicates little difference, although theoretical ultimate moment values 

for the smaller span are about 2 percent higher. 

Table 5.2 and Figure 5.3 show a significant difference between 
experimental and theoretical curvatures at local flange buckling for the 
University of Alberta tests. This may be partly due to the fact that 
curvature in the inelastic range is greatly affected by relatively smal] 
increase in moment, and partly due to the fact that there is a difference 
between the actual stress-strain relationship and that used in the analysis. 
A stress increment of 3 Ksi at a stress level of 60 Ksi for a yield 
stress of 44 Ksi may produce a strain increment of 0.01 based on an in- 
elastic tangent modulus of 300 Ksi. This strain increment of 0.01 cor- 
responds to the strain at the onset of the strain hardening, or to a 
curvature increment 0.00133 for a 12 inch depth section such as a W12 x 36 
section. For example, a beam which buckled locally at a stress of 60 Ksi 
and a curvature of 0.0088, such as Beam 18 in Table 5.2, may show a 15 


percent difference in curvature in spite of only a 5 percent difference in 
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moment. 
As shown in Figure 5.3, curvatures at the onset of local 
buckling decrease fairly rapidly with increasing amounts of longitudinal 


slab reinforcement. 


5.3 Lateral Buckling 
5.3.1 Presentation of Data 


As indicated in Chapter II, lateral buckling of composite beams 
under negative moment may occur for beams where local buckling is pre- 
cluded. Beams tested by Piepgrass were composed of a W12 x 16.5 steel 
section with a 3" x 3/8" coverplate welded on the compression flange and 
a concrete slab with various amounts of longitudinal reinforcement as 
Shown in Figure 5.5, Table 5.4 compares Piepgrass' test results with 
analytical values. Figure 5.6(a) shows analytical values of M/M, for 


various span lengths and amounts of longitudinal reinforcement. The 


effect of varying size of coverplate is shown in Figure 5.6(b). 
DEO se DVESCUSS ION 


The difference between experimental and theoretical ultimate 
moment values are within 5 percent as shown in Table 5.4. This agree- 
ment indicates that it is acceptable to assume a hinge located at the 
end of the stem of the inverted tee section. 

The effect of span length on lateral buckling values is 
shown in Figure 5.6(a) for the section used in Piepgrass' test beams. 


The ratio M/M decreases with increasing span, and the ultimate moment 
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values become less than the plastic moment values in the span range of 
115 to 165 inches with longitudinal slab reinforcement varying from 
0 to 2.4 in? The length of the negative moment region for actual beams 
may vary from 50 to 150 inches. Therefore lateral buckling effects 
must be considered when designing continuous beams. 

Figure 5.6(a) shows the effect of coverplate size on lateral 
buckling. A wider coverplate produces a lower value of M/M,: 

Figures 5.6(a) and (b) and Table 5.4 show that the ratio 


mae, is influenced by the amount of longitudinal slab reinforcement. 


5.4 Moment-Curvature Relationships 


Moment-curvature relationships for various beam sections are 
shown in Figures5.7 to 5.13. Figures 5.7 and 5.8 show the effect of 
residual stress patterns on moment-curvature relationships in a positive 
moment region. Residual stresses at flange tips are assumed 0, 0.3 and 
0.50, Residual stress pattern I, which is a linear stress pattern, 
produces little effect on the moment-curvature relationships. However, 
residual stress pattern II, which is a parabolic pattern, produces con- 
siderable effect on the relationships as shown in Figure 5.8. Test results 
shown in these figures are those obtained by Ferrier and Davison. The 
test results for a plain W12 x 36 section compares favorably with the 
analytical results based on pattern II with tes 0.30, The experimental 
moment-curvature relationships for composite sections do not fit the 
analytical relationships in the plastic range, but they compare favorably 


in the elastic range. The presence of residual stresses decreases the 
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stiffness of the composite sections in positive bending as shown in 
Arguvestie2 anded28 ; 

Figures 5.9 and 5.10 show the effect of residual stress on 
moment-curvature relationships in a negative moment region. These 
figures include results of Davison's tests for comparison. With in- 
creasing amounts of longitudinal reinforcement, the actual stiffness of 
the composite section decreases in comparison with theoretical stiffness. 
This may be due to the effect of slip. According to Davison's fect ; 
the interaction factor, which is discussed in Chapter III, was approximately 
equal to 0.6 when the longitudinal reinforcement reached yield strain. 

From a comparison of theoretical and test results, the interaction factor 
appears to decrease with increasing stress in the longitudinal slab 
reinforcement, since theoretical values based on the interaction factor 
agree with test values in the elastic range, but not in the plastic range. 
Presumably this difference may also relate to the behavior of shear con- 
nectors. 

The effect of residual stresses on the moment-curvature rela- 
tionship in a negative moment region differs from that in a positive moment 
region. Since the stress in the bottom flange is compression in nega- 
tive bending and tension in positive bending, residual stresses have a 
different effect on moment in positive and negative bending. The stiff- 
ness is increased by residual stresses after yielding in a negative moment 
region as shown in Figures 5.9 and 5.10. 

Figures 5.11 to 5.13 show experimental and theoretical 
moment-curvature relationships for beams tested in the present study. 


Theoretical relationships are in good agreement with experimental rela- 
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tionships in positive moment regions. However the agreement is not good 
in negative moment regions, except for Beam CB1 which has the smallest 
amount of longitudinal slab reinforcement. In the elastic range, the 
difference is large with the actual stiffness approximately two times the 
theoretical stiffness. Strain-hardening begins shortly after initial 
yielding of the steel in the negative moment regions. All three beams 
tested ‘attained ultimate moments greater than simple plastic moments 


as shown in Figures 5.11 to 13. 


5.5 Deformation 
5.5.1 Deflection 


Load-deflection relationships obtained in previous investiga- 
tions and in the present study are compared with theoretical relation- 
Ships in Figures 5.14 to 5.19. Figure 5.14 compares theoretical and ex- 
perimental load-deflection relationships for a simple span tested by Yam 


(41) | Yam's analysis includes the effect of slip. The load- 


and Chapman 
deflection relationship based on the writer's proposed analysis dis- 
cussed in Chapter III is also shown. In the elastic region the proposed 
analysis agrees well with Yam's analysis. As ultimate load is approached, 
the proposed method results in deflections larger than those obtained in 
Yam's analysis. This difference is due to the fact that the stress-strain 
relationship for concrete employed by Yam is an idealized elasto-plastic 
relationship whereas that employed in the proposed analysis is based in 


Desayi's approximate formula as given in Equation 2.8. The latter rela- 


tionship has a falling portion after maximum stress, causing reduced 
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stiffness at high loads. Both Yam's analysis and the proposed analysis 
yield values close to experimental results. 

Some complete load-deflection relationships for continuous 
beams have been reported in previous investigations (7) (11) (12) | Test 
results presented by Culver et ai 7) are compared with theoretical values 
in Figure 5.15. Theoretical values are in good agreement with the test 
values. 

Figure 5.16 compares results of tests on continuous beams 
conducted at Lehigh University by Daniels et ai (11) with theoretical 
values. Since the beam span was comparatively long, the effect of shear 
on the deflection is relatively small. However, it still amounts to 17 
per cent of the bending deflection. Theoretical results are in fairly 
good agreement with experimental results. In the simple plastic analysis 
hinges form almost simultaneously at the interior support and load points. 
The load point hinges form first, followed closely by the formation of 
the hinge at the interior support. 

Deflection at the load point at the formation of the last 
hinge is determined by the slope-deflection method. The coefficients 
of the beam stiffness matrix for a composite beam are related to the ratio 
of the stiffness in positive bending to that in negative bending and are 
also related to the ratio of the length of the positive moment region 
to that of the negative moment region. Appendix C includes tables of 
values for stiffness matrix coefficients for the slope-deflection method. 
The load-deflection relationships for beams tested in the present investi- 
gation are compared with theoretical relationships in Figures 5.17 to 


5.19. The effect of shear on deflection is appreciable, since the beams 
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are relatively short. As indicated in the previous section there is 

a relatively large difference in the moment-curvature relationships in 
negative moment regions for Beam CB2 and Beam CB3. However, the analytical 
load-deflection relationships are in good agreement with the test rela- 
tionships. In the analysis the effective width is taken as 40 inches 

as discussed in Chapter IV. 

In the simple plastic analysis the first hinge forms at the 
interior support in all beams. The second hinges at the load points form 
immediately after the first hinge formation for Beams CB1 and CB3. 

For Beam CB2 the first hinge forms at a load of 95 Kips and the second 
hinge forms at a load of 121 Kips. Deflections in the elastic range, 


determined by simple plastic analysis, are small. 
5.5.2. Rotation 


Little load-rotation information for continuous composite beams 
is available in previous research. A certain amount of rotation data was 
obtained in the present study. As discussed in Chapter IV, rotations 
were measured at the ends of beams. Results are shown in Figures 5.20 
to 5.22 together with analytical relationships. The characteristics of 
the load-rotation relationships are similar to those of load-deflection 
relationships. The effect of shear on the end rotation is not as sig- 
nificant as on the center deflection. The ratio of the center deflection 
due to shear to that due to bending is 12k_, E1/AGL? for a simply supported 
beam with span length 2 and a concentrated P at the midspan, and the ratio 
of the end rotation due to shear to that due to bending is 8k, EI/AG2". 
This implies that the effect of shear on end rotation is two-thirds of 


that on the center deflection. 
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5.6 Load-Moment Relationships 


The theoretical load-moment relationship is linear up to 
initial yielding. After yielding begins the moment is redistributed and the 
relationship is no longer linear. Figure 5.23 shows experimental and 
theoretical load-moment relationships for beams tested in the present 
investigation. 

The theoretical load-moment relationships are in good agree- 
ment with test results, with the theoretical values being slightly greater 
than the test values for positive moment in the elastic range. This in- 
dicates that the stiffness of the beam in a positive moment region was 
greater than theoretical stiffness, or the stiffness in the negative 
moment region was smaller than the theoretical stiffness. 


For a two-span continuous beam with a load P at each midspan 


tee M,/2 = PL/4 en 


where M,and M, are bending moments at the load point and the interior 
Support, respectively. The increment in positive moment, AM, there- 
fore, produces an increment in negative moment equal to 2AM. The dif- 
ference between the theoretical and the test load-moment relationships 
is large in the negative moment region. 

The experimental values of the negative moment ve in Beam CB2 
began to increase fairly rapidly at a load of 100 Kips and to approach 
the positive moment M- This implies that the ratio of stiffness in the 


negative moment region to that in the positive moment region increased for 
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loads greater than 100 Kips. This condition can also be seen in the moment- 
curvature relationships for Beam CB2 in Figure 5.12. 

The moment redistribution in Beam CB3 commenced at a load 
of 50 Kips and thereafter negative moment ceased to increase proportion- 
ally with the load due to a decrease in stiffness. This behavior can be 
explained by the moment-curvature relationships for Beam CB3 shown in 
Figure 5.13. The test curve in the figure shows yielding at approximately 
800 in-Kips which is equivalent to the load of 50 Kips as indicated in 
the load-moment relationships in Figure 5.23. Since the stiffness in 
the negative moment region is greater than that in the positive moment 
region after yielding, as shown in Figure 5.13, the negative moment 
increased at a higher rate when the load reached 60 Kips, when the posi- 
tive moment at the load points reached the yield moment. Load-moment 
relationships based on simple plastic theory shows good agreement with 


the experimental relationships. 
5.7 Failure Modes 


Figures 5.24 show the relationships between the failure loads 

and the amounts of longitudinal slab reinforcement for beams tested in 

the present study. Failure modes are local buckling in the negative 
moment region and crushing of concrete in the positive moment region. 

The ultimate concrete strain at failure is determined by Equation 2.9 

as 0.0032 for a concrete strength of 5.5 Ksi. Experimental and theore- 
tical ultimate loads are shown in Table 5.5, Beam CB1 failed due to 
crushing of concrete at a load of 135 Kips, after which local buckling 


occurred at the interior support. The proposed analysis predicts crushing 
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of concrete at a load of 135 Kips, and local buckling at a load of 141 
Kips which of course means that failure is predicted at a load of 

135 Kips. Since the stiffness in the positive moment region decreased 
considerably after crushing of the concrete, local buckling occurred 
at a load lower than the ultimate load. 

Beam CB2 failed in local buckling which initiated at a load 
of 120 Kips and completely formed at 130 Kips and subsequently in crushing 
of concrete at a load of 133 Kips. Figure 5.24(b) shows that analysis 
predicts local buckling at a load of 122 Kips and crushing of concrete 
at a load of 131 Kips. In a continuous beam the ultimate load can be 
greater than that causing local buckling because of redistribution of 
moment. 

Beam CB3 failed in local buckling which began at a load of 
85 Kips and completely formed at a load of 90 Kips, followed by crushing 
of concrete. The analysis predicts crushing of concrete at a load of 
89 Kips and local buckling at a load of 90 Kips as indicated in Figure 
5.24(c). As in the case of Beam CB2, local buckling was not the direct 
cause of failure. 

Theoretical predictions for local buckling and crushing failures 
are in good agreement with experimental results for all beams. Although 
local buckling actually initiated at lower loads than predicted for 
Beams CB2 and CB3, the beams attained the predicted loads based on con- 
crete crushing at an ultimate concrete strain of 0.0032. This may indi- 
cate that the beam stiffness was not significantly reduced by the ini- 
tiation of local buckling. Crushing of concrete causes immediate failure 


and the beam cannot sustain additional load. All three beams tested in 
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the present study attained ultimate loads 8 to 12 percent greater than 
those predicted by simple plastic theory. 

Ultimate loads based on simple plastic theory are not directly 
proportional to the amount of longitudinal slab reinforcement. Values of 
ultimate load do not change significantly for areas of longitudinal slab 


reinforcement in excess of Ao >» Since the plastic negative moment 


yy 
does not increase appreciably once the neutral axis reaches the tension 
flange. The limiting value Ayoy! ony is shown in Figures 5.24. Failure 
loads due to concrete crushing increase proportionally with the simple 
plastic ultimate loads. Ultimate loads based on local buckling are almost 
constant regardless of the amount of reinforcement, although local buckling 
loads for simply-supported beams are almost proportional to the simple 
plastic ultimate loads as discussed in Section 5.2. This is due to the 
fact that an increase in longitudinal slab reinforcement produces an 


increase in negative bending moment capacity because of increased stiff- 


ness. 


5.8 Behavioral Study 


Additional theoretical analyses were conducted for a number 
of composite beams with a steel section varying in depth from 10 and 16 
inches and with flange width-thickness ratios less than 54/voy which is 
a limitation for plastic design of plain steel sections. This limiting 
width-thickness ratio is 8.14 for G40.12 steel with a yield stress of 
44Ksi. The beams were two-span continuous beams with each span equal 


to 20'-0". Slab thicknesses of 4 and 6 inches were chosen. A slab 
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width of 60 inches was used in all beams. The concrete strength was 
3500 psi. The amount of longitudinal reinforcement was varied up to a 
maximum value equal to twice the area of the web of the steel section. 
Steel properties such as yield stress and strain-hardening modulus were 
the same for the steel section and the reinforcement for convenience. 
Details of the beams are shown in Figure 5.25. 

Analytical results related to failure loads and failure modes 
are shown in Figures 5.26 and 5.27. The failure modes are local buckling 
at the interior support and crushing of concrete at the load points. 
Ultimate concrete strains of 0.003, 0.0035 and 0.004 are considered in 
determining ultimate loads based on a crushing failure. The ultimate 
loads are expressed by the ratio to the ultimate simple plastic load for 
the beam with no longitudinal slab reinforcement on the vertical co- 
ordinate and the amounts of longitudinal slab reinforcement are expressed 
by the ratio to the web area of the steel section on the horizontal 
coordinate in the figures. 

The analyses indicate that the flange width-thickness ratio 
significantly effects the ultimate loads based on local flange buckling. 
W12 x 27 and W16 x 36 sections have b/t ratios equal to 8.14 which is 
equal to the upper limit of 54/ Voy. For these two sections the local 
buckling loads are very close to the simple plastic ultimate loads for 
longitudinal slab reinforcement greater than A and the local buckling 
loads are smaller than the simple plastic ultimate loads for amounts of 
longitudinal reinforcement equal to 2A The local buckling loads for 
such stocky sections are significantly greater than the simple plastic 


ultimate loads. 
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Ultimate loads based on concrete crushing increase almost 
proportionally with the increase in the simple plastic ultimate loads. 
Ultimate load based on an ultimate strain of 0.0035 is greater than the 
simple plastic ultimate load by more than 5 percent. Ultimate load 
based on a strain of 0.003, the value employed in the ACI 318-71 Build- 
ing Code, is greater than the simple plastic ultimate load by only 2 to 
3 percent. 

The failure load curves for W12 x 27 in Figure 5.26 are 
quite similar to those in Figure 5.24(b), in spite of different beam 
length, different strength of concrete and different slab width. Results 
of the analysis for 4 inch and 6 inch slab thickness are very similar. 

Ultimate loads for composite beams with W10 x 25 and W14 x 34 
sections are greater than theoretical values based on simple plastic 
theory for an amount of longitudinal reinforcement equal to 2A. The 
flange width-thickness ratios for W10 x 25 and W14 x 34 sections are 
6.5 and 7.4, respectively, or 20 and 10 percent less, respectively, than 


the provision in CSA Standard S16-1969. 
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TABLE 5.1 COMPARISON OF EXPERIMENTAL AND THEORETICAL ULTIMATE 


MOMENT VALUES FOR BEAMS TESTED BY DAVISON‘>) ano Lever ®) 


AREA OF 
BEAM | STEEL {LONGITUDINAL 
NO. |SECTION }|REINFORCEMENT 
(in?) 
18 
im! 
12 |W12x36 
ils 
14 
21 
22 
23. «| W12x31 
24 
25 
3) 
32 
W12x27 
33 
34 


ULTIMATE MOMENT 
in-kips) 
ANALYSIS _ 
BASED ON 


LAY'S PROPOSED 
SHEAR MODULUS 


EXPERIMENT BIJLAARD'S 
PROPOSED SHEAR 


MODULUS 
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TABLE 5.2 CURVATURE AT LOCAL FLANGE BUCKLING MOMENT FOR BEAMS TESTED BY 
pavrson(®) anp LeveR‘®) 


CURVATURE 
AREA OF (radians/inchx10° 
BEAM STEEU LONGITUDINAL ANALYSIS 
NO. SECTION REINFORCEMENT BASED ON 
(in2) EXPERIMENT LAY'S SHEAR 
MODULUS 


W12x36 


W12x31 
W12x27 


8. 
8. 
6. 
4. 
2) 
4. 
3] 
a, 
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TABLE 5.3 COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS FOR BEAMS 


TESTED BY CLIMENHAGA‘!®) 


RATIO OF 


STEEL ICOVER| BEAM | LONGITUDINAL EXPERIMENTAL 


SECTION LENGTH 


STIFFENER TO ANALYTICAL 
ULTIMATE 
MOMENT 


5x3RSJ 
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TABLE 5.4 


STEEL 
SECTION 


W12x16.5 
W12x16.5 
W12x16.5 
W12x16.5 


LATERAL BUCKLING MOMENT FOR BEAMS TESTED BY prepgrass \*) 


RATIO OF 
EXPERIMENTAL 
TO ANALYTICAL 


AREA OF SIMPLE ULTIMATE MOMENT 
LONGITUDINAL PLASTIC (in-kips) 


REINFORCEMENT | MOMENT 
(in2) (in-kips) [EXPERIMENT | ANALYSIS SaGHERCE 
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CHAPTER VI 


SUGGESTIONS FOR ULTIMATE STRENGTH DESIGN 
OF CONTINUOUS COMPOSITE BEAMS 


6.1 Introduction 


At present no design standards specifically cover ultimate 
strength design of continuous composite beams, although theoretical 
and experimental investigations provide substantial information on ul- 
timate strength behavior of such members. British Standard CP 117 Part 
1 includes provisions for ultimate strength design of simply-supported 
composite beams. However, additional requirements may be necessary for 
ultimate strength design of continuous beams. Such further requirements 


as discussed herein, are mainly based on the present investigation. 


6.2 Present Specifications 


In British Standard CP 117 Part 1, which specifies require- 
ments for simply-supported beams, the specified load factor is 1.75 
for both dead and live loads. The stress in the steel section at ulti- 
mate is the specified yield stress, whether in tension or compression. 
The concrete compression strength is 4/9 of the specified concrete cube 
strength. This fraction is derivéd from the assumption that the strength 


of the concrete in the slab is 2/3 of the cube strength and the pro- 
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vision of a higher factor against concrete crushing than the factor against 
yielding of the steel. This latter provision results in assuming com- 
pression strength as 2/3 of the slab strength. Furthermore the standard 
requires that in no case should elastic stress under working loads in 

the steel section exceed 0.9 times the specified yield stress and in 

the concrete exceed one third of the concrete strength. Deflections are 
calculated on the basis of a fully composite section, using a modular 
ratio of 15 for live loads and 30 for dead loads. 

Reference might be made to ACI 318-71 for concrete ultimate 
strength design requirements in a positive moment region and CSA Stan- 
dard S16-1969 for steel design requirements in a negative moment region. 
Provisions of ACI 318-71 for flexural members permit the use of an 
equivalent rectangular concrete stress distribution in which a concrete 
stress of 0.85 fies is assumed uniformly distributed over an equivalent 
depth of compression zone. Ultimate moment capacity is modified by 
introducing a capacity reduction factor equal to 0.9. CSA Standard S16- 
1969 requires a maximum flange width-thickness ratio for plastically 
designed steel sections. It also provides for the design of shear 


connectors based on ultimate strength considerations. 


6.3 Moment Capacity 


(8) (20) (24) indicate that ultimate moment capa- 


Test results 
city in a positive moment region can be satisfactorily evaluated on the 
basis of an idealized stress distribution in which the steel stress is 
equal to the yield value and the stresses in the concrete slab are repre- 


sented by an equivalent rectangular stress block with a value of stress 
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Test results (2) (6) (10) indicate that ultimate moment capa- 
city in a negative moment region can be satisfactorily evaluated on the 
bases of an idealized stress distribution in which stresses in the steel 
section and in the longitudinal slab reinforcement are equal to their 
respective yield values. However, it is necessary to check for the 


possibility of premature local flange buckling which may have the ef- 


fect of reducing the moment capacity. 


6.4 Flange Width-Thickness Ratio 


For plastic design CSA Standard S16-1969 requires a maxi- 
mum flange width-thickness ratio equal to 54/vo, for compression flanges 
of plain steel beams. However, local flange buckling of composite beams 
is related to the area of longitudinal slab reinforcement as discussed 
in Section 5.8. Therefore it is suggested that the flange width-thickness 
ratio be a function of the longitudinal slab reinforcement area. For 
composite beams with an area of longitudinal slab reinforcement less 
than the web area of the steel section, a maximum flange width-thickness 
ratio equal to S4/voy is suggested. If the area of longitudinal slab 
reinforcement is greater than the web area but less than twice that 


area, the flange width-thickness ratio should be reduced to 49/70. 


6.5 Shear Connectors in a Negative Moment Region 


CSA Standard S16-1969 requires that shear connectors in a 


negative moment region resist a horizontal shear force equal to Aon 


which represents the force in the longitudinal slab reinforcement at 


laste ont ni eoeeaste Ast notudiaser saste bast Isebi aay wl esa | 
ited ot Isyp9 546 snomestont on def2 [entbudtenol ang nt a rotsas | : 
ots yot Anes oF YiseesosR at tt .revowol . 29uT sv bfety evitanqess i 

-to ot sved yom siofdw patfaoud aprstt Tsao! ows 6moq 70 stttatezo 4 
.vitosqso tnemom only enroube7 to sost 


ae ei 


ofgsf 2eendotiT-dibtW spnsli 4.3 
spree 


-ixem 6 esvtupey C8@l-Ol2 bisbna¥2 Aed netesb ofseslg 107”! Whe 
2epnait notezarqnoo ‘tot yore ot [sups offer azaniotds-ridbtw pret? mum 
omesd a¥izoqno> to ent hloud sedslt feool .revewoH .emsed: fast2 nistq to 

bagevoetb 26 snamaovetatey defe Tantbustedof to sov6 aft os basarennet 

e2mtaidd-Atbiw spnelt oad tart bed2oppue ef +1 evoTsredT .8.2 tiofsos2 nt 
4oi . #048 EnamsovoTatey dele Tenfbutipnol sfd to notsanat 6 od otis 

cool tramaovotntey dsfe fantbuttpaol to se16 ne nstw emsed 9s t2oqmoo 
azondotrt-tbiw spastt muntxem @ .noftose Peete edt to ses dow ont not 
dsl2 Jantbusipno! to seve sad tI .betespoue 2? Poa ot [sups offer 

seit sotwi neds east dud sors dow ond nent r9tsexp 2t jnomaovotntet 


gow ven os baoubsy od bfuode ofde: eesndotdd-dthtw sonelt ‘odd .5978 


notpal tnamoM evitspaK 6 nt evotsennod ysod2 €.3— 
; ; - ah. ><. rr 


alii _¥ 


6 Ni evotosnroD ‘sere sods eortupst Qd0r-ar2 byebast2 ACD 8 
yeah 8 feups ato? “sere Tsdnostod s Jereey motpon tnemom svttsgen 


ts dnamaavotnioy dele Ientbustgno! ent ni eovot ods Come a 


188 


yield conditions. Based on results of pushout tests Van Dalen et ai Ika) 


proposed a design load per connector in a negative moment region equal 

to 80 percent of that in a positive moment region. However, results of 
beam tests conducted by Datacon es and lever and results of the pre- 
sent beam tests do not indicate that a reduction in capacity is indicated. 
It is therefore suggested that the design of shear connectors in a nega- 
tive moment region be based on the same connector strength as for a 
positive moment region, provided that the longitudinal slab reinforce- 


ment is adequately anchored in a positive moment region. 


6.6 Deflection 


At working loads, the deflection of composite beams may be 
based on elastic theory. The bending deformation at midspan may be 


expressed as 


for a concentrated load P at midspan. The coefficients k and ko are 
functions of the ratio of stiffness in the negative moment region to 
that in the positive moment region as discussed in Appendix C. Values 


of these coefficients are given in Figure C.2(b). For a simply-supported 
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beam KY and Ko are equal to 5 and 4, respectively, and for fixed end 
beams with equal stiffness in positive and negative moment regions 
kK, and Ko are each equal to 1.0. 

Shear deflection of composite beams may be significant as 
discussed in Section 3.5. Since shear is primarily resisted by the web 


(67) 


of the steel section only , the shear deflection may be expressed as 


for a concentrated load at the midspan. 
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CHAPTER VII 


SUMMARY AND CONCLUSIONS 


7.1 Summary 


Behavior of continuous composite beams has been investigated 
theoretically and experimentally. Analyses for local flange buckling 
and lateral buckling in a negative moment region are proposed and com- 
pared with experimental results. Deformations including the effects of 
shear and slip were studied in the elastic and inelastic region. Three 
two-span continuous composite beams were tested in order to provide addi- 
tional information on deflections, moment redistribution and ultimate 
loads. The test beams varied in terms of steel section size and amount 
of longitudinal slab reinforcement in the negative moment region. Behavi- 
oral studies were conducted in order to obtain failure loads and modes 
of failure for continuous beams in which steel section size, concrete slab 
thickness and amount of longitudindal slab reinforcement were varied. 
Based on theoretical and experimental results, requirements are proposed 


for ultimate strength design of continuous beams. 


7.2 Conclusions 


The major conclusions resulting from the present investigation 


ger 


le The proposed analysis for local flange buckling in a negative 
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moment region gives good agreement with test results. The ratio of 
local flange buckling moment to simple plastic moment is decreased signi- 
ficantly with increase in the flange width-thickness ratio and is slightly 
affected by the amount of longitudinal slab reinforcement and the span 
length. Curvature at the local flange buckling moment is decreased sig- 
nificantly with increase in the amount of longitudinal slab reinforce- 
ment. 

the The proposed lateral buckling analysis based on thin-walled 
beam theory for composite beams in a negative moment region agrees with 
results of beam tests. The analysis indicates that the ratio of lateral 
buckling moment to simple plastic moment is decreased significantly with 
increase in span length and is slightly affected by the amount of longi- 
tudinal slab reinforcement and slightly affected by the size of coverplate 
on the compression flange. 

3s The proposed analysis provides satisfactory predictions of 
ultimate load and failure modes, i.e., crushing of concrete in a positive 
moment region or local flange buckling in a negative moment region. 

4, Failure modes are significantly affected by the amount of 
longitudinal slab reinforcement in the negative moment region. 

ie). The proposed analysis for deflections, which includes the 
effect of shear and slip, predicts satisfactorily the actual deflections 
obtained from the tests. The shear deformation was significant in the 
beams tested. 

6 The proposed analysis for evaluation of moment-curvature 
relationships based on Newmark's integration is in good agreement with 


experimental values. 
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Whe The effect of residual stress on the moment-curvature relation- 
ships in a negative moment region differs from that in a positive moment 
region. 

8. The assumption that slip strain is constant throughout the 
shear span is a satisfactory approximation in the evaluation of the effect 
of slip on deflection. 

9. A maximum compression flange width-thickness ratio equal to 
54/voy is required when the amount of longitudinal slab reinforcement is 
less than the web area of the steel section in order to prevent local 
flange buckling before a mechanism forms. This maximum ratio is reduced 
to 49/Voy for an amount of longitudinal slab reinforcement greater than 
the web area but less than twice the web area. 

10. Moment redistribution in continuous beams is affected signi- 
ficantly by the inelastic stiffness in the positive and negative moment 
regions. 

1 Shear connectors in a negative moment region are capable 
of developing the full yield value of the longitudinal reinforcement 
if sufficient development length of the longitudinal slab reinforcement is 


provided into the positive moment regions. 
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APPENDIX A 
STABILITY EQUATIONS FOR THIN-WALLED BEAMS 


A.1 Kinematics of Deformation 


If displacements —€ and n in x and y directions, respectively, 


and angle of twist 6 are applied to a thin-walled beam section at A(a,»a,)) 


as shown in Figure A.1, the displacements EB and np of an arbitrary point 
B(b, sby) are 


= - (by - a )o A.1(a) 


Np = 7 a )e A.1(b) 


where the sign convention for 6 is based on the right hand rule. The 


displacements v_, We in the tangential and normal directions to the 


S 
cross section at S are 


< 
iT] 


E.cosa + n,sina A.2(a) 


Ww. = -E,Sina + n.cosa A.2(b) 


where Es and ng are displacement at S in the x and y directions, and «a 
is the angle between the tangent at S and the Ox axis as shown in Figure 


A.1. 
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Substituting Equations A.1 with respect to S into Equations 


A.2, the displacements at S in the tangential and normal directions become 


V, = §coSa + nsina + or, A.3(a) 

Wes Seas Nadine COSaet er, A.3(b) 
where ay fe (x - a, ) Sina - (y - ay) cosa 

Cee (x - ay) cosa + (y - ay) sina 


and are shown in Figure A.1. 
Noting that the shear strain ee is assumed to be zero for 


a thin-walled beam section in Figure A.2(49)(50) , 


ou OV 
as SSS, oa ——_ = . 
Usa BS Ye 0 A.4 


the displacement u is obtained by integrating Equation A.4 


u=c- [as A.5 


where < is the initial displacement in z direction. 
By differentiating v with respect to z in Equation A.3(a) and 


multiplying by ds, Equation A.5 may be expressed by 


u=c-e'x -n'y - 0'w A.6 
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since ds cosa = dx 
ds sina = dy A.7 
ry ds = dw 


The quantity w is independent of x and y coordinates, and is 


called the warping function. The strain EF in the z direction is obtained 


by differentiating Equation A.6 with respect to z. 
ce = a ch ae ms n'y = By A.8 


A.2 Normal Stress-Strain Relationship 


From Hooke's law and E. equal to zero, the stress-strain 


relationships with respect to z and s directions may be given by 


From these two equations, the relationship between the normal stress o, 


and normal strain e, may be obtained by 


where E = ae 
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E is defined as the modified elastic modulus. Since v2 is small, it 


may be neglected and the elastic modulus E used instead of E. 


the normal stress becomes 


Spee aX a ve On) 


A.3 Equilibrium Equations for a Thin-Walled Beam 


Then 


A.9 


For a small element shown in Figure A.3 equilibrium in the 


longitudinal direction requires that 


Equilibrium of a strip of thin-walled beam between z and 


(z + dz) as shown in Figure A.4 may be expressed by the following equa- 


tions 


d0 
= Z a . 
LZ S80: | gt beds + (T) Tp a q,) dz = 0 
aT 
=X = 0, | se~ tdz cosa ds + q,dz = 0 
S 


OT 
rY = 0, | = tdz sina ds + qydz = 0 
S 


JZ 


+ Mie dz a+ m,dz = 0 


A.11(a) 


A.11(b) 


A.11(c) 


OT oS 
=M, a= 03 tdz E ~ a.) Sina - (y - a) cosa|ds 
S 


A.11(d) 
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where Meg is Saint Venant's torsional moment defined as 


Mig = GJe 


in which J is St. Venant's torsional constant. 
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Substituting Equations A.7 into Equations A.11 results in 


A.13(a) 


A.13(b) 


Ay i3{c) 


A.13(d) 


The first term in Equation A.13(d) contains the stress re- 


Sultant due to the warping torsional moment. Differentiating equations 


in Equation A.13(b), (c) and (d) part by part and noting that shear flow 


ae t at the boundary equals the external force, that is, 


yields Equations A.14 as follows: 
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ele} 
Z — 
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The functions 30/22 and at. 5/ 98 in Equations A.14 are ob- 


Equations A.9 and A.10 as 


00 
a= FE Ga! su aly ~ n'y _ 6"w)! A.15(a) 
Z 
oT 00 
SRA = OL a *, Ihe " = " a4 " \ 
ee 2 oF (c' - E"x - n"y - 6"w) A.15(b) 


Substituting Equations A.15 in Equations A.14, the following 


equations result: 
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Ee ee ut " AV es ! - 
( 5x5 Lye 4 In : Bue ) dy i TY ThYR 


AN G(s) 

EN fe el ec pe nee = 

Mm, oF Tyo ~ TpWp A.16(d) 

where 
A = | dA 
A 

Sy = | ydA, Sy = i xdA, S = |. wdA 
ly = | xydA, YY = | XwdA, Sy = | wydA A.17 


Equations A.16 may be simplified by choosing a particular 
coordinate system, i.e., the principal coordinates. Some of the values 
in Equation A.17 become zero for principal coordinates and Equations 


A.17 may be expressed by 


E (Ac')' = Gh Tee Tp A.18(a) 


SNP au, Pky SU A.18(b) 
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E (Iyn ee Aye TY > Tpye A.18(c) 


E (16")" - G (Je")" =M, + Tio - Trup A.18(d) 


By introducing Equation A.9 the stress resultants in the 
form of a normal force She = | oda, bending moments MY = [opydn and 
A 


My = [o,xaa, and a bimoment By = fo Jud are defined as 


N#= EAc A.19(a) 
Me aerl ies A.19(b) 
the aak a ae ce 
Beg etnOn A.19(d) 


The normal stress Oo. may be expressed in terms of stress 


resultants by substituting Equations A.18 into A.9 to give 


= 
EX 
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SEE ked Teeth A.20 
A pe ell, 
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A.4 Stability Equations for a Thin-Walled Beam 


In order to investigate torsional buckling of a thin-walled 
beam subjected to an axial force NS and bending moments iy and Ms equi- 
librium after a small displacement must be stablished. The normal stress 
s due to an axial force N, and bending moments MY and My is given by 


Equation A.20 as 


= 


NM 
Gena rey egneay A.21 
i i I 
y X 


Referring to Figure A.5, the components of the internal 


stress oie in x and y directions are 


and 


on oo, us aon. 
a + —2 fea, 
GAA a5 7 the) INAS pores or resp Ee 


respectively. Neglecting smaller quantities, the components can be ex- 
pressed by (0,82) 'dzaa and (a,ng)dzaA, respectively. The total trans- 
verse forces per unit length qy and qy due to the components can be 


obtained as 
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sete I, (cé2)' dA A.22(a) 


qy = f (cer uadn A.22(b) 


given by 
aes |. { (o,f) iar Bh ale eo aaly - ay)} dA A.23 


Substituting Equations A.1 into Equations A.22 and A.23, 


and integrating across the cross section, results in 
she {N,(e + 20) SEP A.24(a) 
dy = {Nn - a6 )} + Hous A.24(b) 


A.24(c) 


where M, is a new stress resultant defined by 


H Joa? Eat yee ay)?} A.25 


Stability equations for thin-walled beams result from substituting 
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at \p bne 2 hk ang pe 
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Equations A.24 into Equations A.18 as 
ESA eave = {Ce + aye)} - (M,o')' + Tix, - Toxp A.26(a) 
E (1yn")" = {ti,(n' - aye')} + (Wyo")! + Ty, - Teyp A.26(b) 


E (Lee) - G(Je')' = {N, (-aya" + aye") - a OP s 


G clietadl es y EIB A.26(c) 
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S(x,y) 


FIGURE A.1 DISPLACEMENTS FOR A THIN-WALLED BEAM 
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FIGURE A.2 CO-ORDINATE SYSTEM 


FIGURE A.3 


STRESSES IN LONGITUDINAL DIRECTION 
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FIGURE A.4 


STRESSES ON AN ELEMENT OF A THIN-WALLED BEAM 
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FIGURE A.5 DISPLACEMENT OF A STRIP 
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APPENDIX B 
APPROXIMATE METHOD FOR DETERMINING THE EFFECT OF SLIP 


The effect of slip on the stiffness of a composite beam 
may be calculated on the assumption that slip strain is constant along 
the shear span and is defined as the slip divided by the shear span. The 
stress-strain relationships for concrete and steel are assumed linearly 
elastic, the tensile strength of concrete is assumed to be zero and the 
force-slip relationship is assumed to be linear. 

There are two specific cases to be considered in the evalu- 
ation of the stiffness as in the evaluation without slippage; namely, 
the neutral axis is in the concrete slab (Case I) and the neutral axis 


is in the steel section (Case II). 


Case I: Neutral Axis in Concrete Slab 
Based on conditions shown in Figure B.1(a), the equilibrium 


equation for the beam is 


| odA = | odA + | oda 
A ns A. 


b 
~$E, 5 V9 + 6 EAL (F+ to - yy - eg/o) = 0 Bel 


where Eq is slip strain and can be expressed as 
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= A/L, Be 


and L. is the shear span. The relationship between concrete force C 


and slip A. is expressed by 
Ee Ga Bea 


The shear force on the shear connectors in the shear span 
is gE by4/2. Slip strain is then obtained from Equations B.2 and B.3 
as 
2 
ee bE DG 
d 2kiv 


SS 


B.4 


By introducing Equation B.4 into Equation B.1, the following 


equation related to Ya is obtained. 


2 a = 'y2)\) = 
Ae ee RENE Pa eae ep ete ESD 
EA nA 
Faee Se Sess 
where p vg Se 
c 

and Ko ed. - ts 

S 2b.L 2nk.L. 


The solution of Equation B.5 gives the location of the neutral 


axis as 
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Yq = TKD (/T + (a + 2b) + Kip)/p - 1) B.6 


For a beam with complete interaction 
‘p= Cre CR aay 20) 


The bending moment is defined by 


=< 
i} 


& | oydA + oydA 
A A 


where Tag is an equivalent moment of inertia for the section and is 


defined as 
= 2 U 
I Ie 3 Ace. tf Is B.8 


in which I. and le are the moment of inertia of the steel section and 
concrete slab, respectively. Le is evaluated by by4/3 and e. is expressed 
:. > 2 
by d/2 + t. Yq Ebyg/2kL.. 
Case II: Neutral Axis in Steel Section 


Based on Figure B.1(b) the equilibrium equation is 


| odA = | odA + | odA 
A A. Ae 
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oE 
G 
aa + 2y4) + gEA. (d/2 Sab eq/o) = 0 B.9 


Introducing the slip strain AiG eb (te * 24) A Jel ok. 


the neutral axis location defined by Yq can be obtained as 


ia = i B.10 
2 mei (1 + kA, 


For a beam with complete interaction 


tA /M - dA. 


py 2A Ty + AG) B. 1] 
The bending moment Myis expressed as 
uh, = Joyan 
= El ag? Bele 
where Tag is an equivalent moment of inertia expressed by 
lag Sela Ae? at We - Nae Baas 
where Et. a 2vqg)Ac 


e =d/2-y, - —S——*— 
s d 20K 


26 be2ee 1qxs ef ft aiach nines 


Abyo| 7 a 


$1.8 


218 


fo ele YG 


In the evaluation of the beam stiffness by the above proce- 
dure, the value of the constant kK. is required. This value may be ob- 


tained from a force-slip relationship obtained from push-out tests. 
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(b) NEUTRAL AXIS IN STEEL SECTION 


FIGURE B.1 EFFECT OF SLIP ON STRESS DISTRIBUTION 
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APPENDIX C 
APPLICATION OF SLOPE-DEFLECTION METHOD TO COMPOSITE BEAMS 


C.1 Deflection at Ultimate Moment for Simple Plastic Theory 


For a segment between sections at which plastic hinges are 
located elastic continuity exists up to the load at which the last hinge 
forms. Therefore slope-deflection equations may be used to determine 
deflections. The slope-deflection equation is given for a beam with 


constant rigidity by 


Ste See eree ed tire Meme 
Snape | AB ty gi a3EL hin | a 
where the sign convention is defined in Figure C.1(a). 

Since stiffness is not constant for a composite beam subjected 
to positive and negative moments, equation C.1 can not be directly em- 


ployed. However, the case may be handled by means of equations similar 


to Equation C.1. 


A I 
ar Q te ae (k Map “i k Mga) Gr2\a) 


8 = 6, 
AB B 
; p 


Sess A 
Sg et 4 + ET, Ck Map + k ga) C.2(b) 


where EI_ is the rigidity of a positive moment region. The quantities 
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Bap and Ben are rotations at ends A and B, as shown in Figure C.1(b), 
due to similar loading of a simply supported beam. The constants kj; 


Kyo, Ko, and ky, which are related to the flexibility co-efficients, 


{1 - 0'(1 - 63)} 


Kio = Ko, = - ae - 9'(3a% - 203) 


are defined by 


Ki 


TT 
ole 


1 ; 
oe tee) 
and o' =1l- 


where ao and g are illustrated in Figure C.1(c) and p is the ratio of 
stiffness in a negative moment region to that in a positive moment region. 
The co-efficients k,;, ky. and ky» are tabulated in Table C.1 to C.3. 

For constant stiffness the rotations Bap and OBn for a uni- 


formly distributed load are given as 


3 
=— 0s = ger Re, 


np ~~°pa = 24ET 


For composite beams with different stiffness in positive and negative mo- 


ment regions, these rotations become 
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where S 1 - p'a2(6 - Ba + 302) 
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1 - p'a3(4 - 3a) 


Co-efficients Sy and So are tabulated for values of a and 0 
in Tables C.4 and C.5. The deflection at ultimate moment can be evalu- 
ated from Equations C.2, since the moments are equal to the plastic mo- 


ment values. 


C.2 Slope-Deflection Equations for Composite Beams 


The equations in Section C.1 are applicable between plastic 
hinges, therefore they are based on one positive and one negative moment 
region as shown in Figure C.1. However an interior span consists of one 
positive and two negative moment regions as shown in Figure C.2. For 


such cases the flexibility co-efficients Si % ko are give by 
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where the co-efficients O12 4,5 By» Bos P, and p, are shown in Figure C.2. 


2 
The values of S, and S, in Equations C.4 are given by 


at (6 - 8a, + 302) - p52 (4 - 30, ) 


: a: (4 - 3a, ) - po 


Solving Equations C.2, the end moments Map and Mon may be 


expressed by 
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and Mir and Mor are fixed-end moments at A and B, and are expressed by 


Map = DB Mye 


BAw °>2- BE 


where Mar and MBF are the fixed end moments for a beam with constant 
stiffness. The values of b, and b, for a uniformly distributed load are 


given by 


The midspan deflection for a uniformly loaded beam with span 2 


and fixed ends may be expressed by 
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For a concentrated load at midspan, the midspan deflection becomes 


A = [kK A ae 
c c 192EI 
p 
where Ko = le 
oy C, 


The values of kK and kK. are shown in Figure C.3(a). For a uniform 


load the bending moment Meg at the ends is given by 


Deane ee 
Mee LL 1 
b CuK 
Dek} 1 Ww 
where Mm, = _ at 9 


where 1. 6. ee 


The values of My and m. are shown in Figure C.3(b) as a function of the 
ratio of stiffness in a negative and positive moment region. The moment 
at midspan is obtained from the equilibrium equation for a statically 


determinate beam, if the bending moment at the ends is known. 
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0.20 
0.25 
0.30 
0.35 
0.40 
0.45 
0.50 
0.55 
0.60 
0.65 
0.70 
0.75 
0.80 
0.85 
0.90 
0.95 
1.00 


0.20 
0.25 
0.30 
9.35 
0.40 
0.45 
0.50 
0.55 
0.60 
0.65 
0.70 
0.75 
0.80 
0.85 
0.90 
0-95 
1.00 


TABLE C.1 
0.15 0.20 
0.8478 0.9840 
0.7192 0.6213 
0.6335 0.7129 
0.5722 0.635% 
0.5263 0.5773 
0.4905 0.5321 
0.4620 0.4960 
0.4386 0.4664 
0.4191 O.4478 
0.4026 0.4209 
0.3885 0.4030 
0.3762 0.3876 
0.3655 0.3750 
0.35€0 0.3620 
0.3476 06.3574 
0.3401 0.3419 
0.3333 0.3333 
0.60 0.65 
71.5813 17.6055 
71.2693 1.2905 
1.0613 1.0778 
0.9128 0.9258 
0.8013 0.8118 
0.7187 €.7233 
0.6453 0.6524 
0.5886 0.5944 
0.5413 0.5460 
0.50713 0.5051 
0.4670 0.47071 
0.4373 0.4397 
0.8713 0.4731 
0.3884 0.3896 
Q.3€60 0.5686 
0.3458 0.3501 
0.3333 0.3333 


VALUES OF CO-EFFICIENT k,, IN EQUATION C.2 


0.25 


1.1042 
0.9775 
0.7830 
0.6912 
0.6224 
0.5689 
0.5260 
0.4910 
0.4618 
0.4371 
0.4159 
0.3976 
0.3815 
0.3673 
0.3547 
8.3435 
0.3333 


0.70 


1.6307 
1.3063 
1.0901 
0.9357 
0.81798 
0.7297 
0.6577 
0.5987 
08.5496 
0.5060 
6.4723 
0.4u14 
0.4744 
0.3906 
0. 2694 
0.3504 
d. 3333 


0.75 


1.6458 
1.3177 
1. 0990 
6.9427 
0.8255 
0.7344 
0.6615 
0.6018 
0.5521 
0.5700 
0.4740 
0.4427 
0.4754 
8.39742 
0. 3698 
0.3506 
0. 3333 


1.3005 
1.@£37 
0. 8975 
0.7824 
0.6960 
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0.3480 
0. 3333 


0.50 


1.6653 
41.3323 
1.171103 
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0.4583 
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Q285 0.3335 
090 053335 
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0.20 0.6213 
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VALUES OF CO-EFFICIENT k,,, 
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0633 5m 
0.3348 
0.3345 
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0.6080 
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0.3476 
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0.3369 
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0.5893 
0.5419 
0.5040 
0.4730 
0.4471 
0.4252 
0.4065 
0.3902 
0.3760 
0.3635 
0.3523 
0.3423 
O233533 


IN EQUATION C.2 


0.85 


Leso2]e 
0.9475 
0.8110 
0.7135 
0.6404 
0.3835 
0.5380 
0.5008 
0.4698 
0.4436 
0.4211 
0.4016 
0.3845 
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0.3561 
0.3441 
0.3333 


0.4548 
0.4245 
0.4042 
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O36 37 
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0.3536 
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0.3464 
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0.3409 
053387, 
0.3367 
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1.3053 
1.0623 
0.9003 
0.7846 
0.6978 
0.6303 
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Oo 322 
0.4953 
0.4642 
0.4375 
0.4143 
0.3941 
0. 3762 
0.3603 
0.3461 
Ma S3s)3) 


0.5000 
0.4583 
0.4306 
0.4107 
0.3958 
0.3843 
0.3750 
0.3674 
0.3611 
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0.3512 
0.3472 
0.3438 
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0.5239 
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0.4558 
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0.4048 
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0.3651 
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0.5552 
0.4997 
0.4627 
0.4363 
0.4165 
0.4014 
0.3888 
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0.3703 
O23632 
0.3571 
0.3518 
0.3472 
0.3431 
0.3395 
0.3363 
OF3338 


1.00 


1.6667 
173233 
1.11171 
0.9524 
0.8333 
0.7407 
0.6667 
0.6061 
0.5556 
0.5128 
0.4762 
0.4444 
0.4167 
Or 3322 
0.3704 
0.35109 
033333 
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0.40 1.0180 
0.45 1.0146 
0.50 1.0120 
0.55 1.0098 
0.60 1.0080 
0.65 1.0065 
0.70 1.0051 
0.75 1.0040 
0.80 1.0030 
0.85 1.0021 
0.90 1.0013 
0.95 1.0006 
1.00 1.0000 
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0.20 2.9008 
0.25 2.4256 
0.30 2.1088 
0.35 1.8825 
0.40 1.7128 
0.45 1.5808 
0.50 1.4752 
0.55 1.3888 
6.60 1.3168 
0.69 1.2559 
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0.75 1.1584 
0.80 1.1188 
0.85 1.0839 
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VALUES OF CO-EFFICIENT s, IN EQUATION C.4 
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1.0408 
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(a) SIGN CONVENTION 
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(b) END ROTATION 
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APPENDIX D 


COMPUTER PROGRAM 


D.1 Introduction 


The following steps are required in the computations for 
deformations and local flange buckling for a given load: 
Ls Compute moment-curvature relationships in positive and 
negative moment regions. Moment-curvature relationships in positive 
moment regions include the effect of slip. 
ras Compute bending moments and bending deflection by employing 
the finite difference method. Rotation and shear are computed from de- 
flection and moment values respectively. Shear deformation is computed 
by the virtual work method. In the inelastic region the bending deforma- 


tion is obtained by an iterative procedure described in Section 3.9. 


Sh Check local flange buckling in the negative moment region. 
4, Compute the maximum concrete strain in the positive moment 
regions 


The computer program was written in Fortran IV and computa- 
tions were carried out on the IBM 360/67 computer at the University of 
Alberta Computing Center. The flow chart in Figure D.1 outlines the se- 
quence of the computations required for the analysis of the continuous 


beams. The program consists of 24 subroutines and 3 functions. 
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D.2 Description of Subroutines and Functions 


READA 
DEBU 


MC 


READB 


DERM 


LOCAL 


MX 
NAPB 


FORCE 


reads and writes the beam dimensions 

computes and writes elastic properties and consists of sub- 
routines NUE, MISB, MISC, SCF and QXx2. 

computes the moment-curvature relationships in positive and 
in negative moment regions. The outline of MC is given in 
FigurcesD eZ. 

reads and writes the beam length, number of spans, locations 
of intermediate supports, loads and load locations. 
computes bending moment, shear, deflection and rotation in- 
cluding effects of shear and slip. The outline of the com- 
putations is shown in Figure D.2. 

computes local flange buckling values. If the first eigen- 
value is greater than 1.0, buckling may occur. LOCAL con- 
sists of subroutines YAMA, SHUKI, TOKYO and PORT. SHUKI 
and TOKYO compute eigen values and were developed by M. Suko. 
YAMA is the memory for the moment-curvature relationship 
and location of neutral axis. PORT computes inelastic tan- 
gent modulus and inelastic shear modulus. 

computes maximum slab strain in each span. 

computes moment for given curvature or curvature for given 
moment. The outline of NAPB is illustrated in Figure D.3. 
computes normal force and moment for a rectangular element 
for a given material (i.e., steel or concrete), neutral 
axis location and curvature. FORCE includes the functions 
STRES and FOC which are the stress-strain relationships for 


steel and concrete, respectively. 
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computes slip deformation of shear connectors. 

computes redundant forces. REACT includes FDM and SMA. 
computes additional external moment for inelastic regions. 
computes shear deformation and includes SHERG 

computes shear energy in a segment. 

solves equations by Gauss elimination. 

writes computational results of forces and deformations. 
computes deflection, rotation, shear and moments by means 


of the finite difference method. 


D.3 Input Data 


Input data required for this analysis is listed in the 


following, and input format is illustrated in Figure D.4. Explanation 


of the symbols for input data are listed as below. 


B 
- 


SiGY 


Flange width for steel section 

Flange thickness for steel section 

Beam depth for steel section 

Web thickness 

Yield stress for steel section 

Ratio of elastic modulus to strain hardening modulus for 
steel beam 

Ratio of the strain at strain hardening to yield strain 
Ultimate strength for steel section 

Number of segments of flange in thickness direction 
Number of segments of web in depth direction 

Number of segments of flange in width direction 


Area of longitudinal slab reinforcement in negative-moment 
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Yield stress for longitudinal reinforcement 

Distance from centroid of longitudinal reinforcement to 
top of steel section 

H for reinforcement 

RSN for reinforcement 

Interaction factor of the strain in longitudinal reinforce- 
ment 

Concrete slab width 

Slab thickness 

Concrete strength 

Ratio of elastic modulus for steel to that for concrete 
Number of segments of slab in thickness direction 

Beam length 

Number of segments in length direction 

Number of intermediate supports 

Location of intermediate support 

Number of load stages 

Location of concentrated load 

Concentrated load 


Uniformly distributed load 
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START 


READ DIMENSIONS OF COMPOSITE BEAM AND 
PROPERTIES OF STEEL AND CONCRETE (READ A) 


COMPUTE SECTION PROPERTIES (DEBU, MP) 
AND MOMENT CURVATURE RELATIONSHIPS (MC) 


READ BEAM LENGTH AND SPAN, LOAD, LOAD 
LOCATION (READ B) 


COMPUTE DEFORMATION INCLUDING SHEAR 
DEFORMATION (DERM) 


COMPUTE LOCAL BUCKLING VALUES AND 
MAXIMUM SLAB STRAIN (LOCAL, MX) 


PRINT OUT COMPUTATIONAL 
RESULTS (CPRINT) 


T 
LAST LOAD STOP 


( ): SUBROUTINE 


FIGURE D.1 MAIN PROGRAM OUTLINE 
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phot ss OR 


(XM JAD0.) MIART2 GALE WUMIKAM 


oe! = 


2 =e F874 Se 


ENTER 


NUMBER OF 
SPAN, 
NM 


NM > 1 


COMPUTE REDUNDANT, X; (REACT) 


REDUNDANT j= Xj.1#(Xj-Xj-1)/30 
CC =i X;-Xj- 71 


COMPUTE BENDING MOMENT AND DEFORMATION 
INCLUDING REDUNDANT (FDM) 


COMPUTE ADDITIONAL EXTERNAL 
MOMENT (MMI) 


COMPUTE SHEAR DEFORMATION 
(DS) 


RETURN 


( ): SUBROUTINE 


FIGURE D.2 OUTLINE OF DEFORMATION PROGRAM (DERM) 
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ENTER 
SET NEUTRAL AXIS 


COMPUTE Ns=fa,odA, Ms=J,,aydA 
FOR STEEL SECTION (FORCE) 


CURVATURE 
20 
T 


COMPUTE SLIP STRAIN 
(SLIP) 


COMPUTE Nr=Jf, odA,Mr={,,odA 
FOR SLAB REINFORCEMENT 
N =Ns+#Nr, M=Ms +Mr 


COMPUTE Ne=J,,odA, Mc=f,.oydA 


FOR CONCRETE SLAB (FORCE) 
N=Ns+#+Ne, M=Ms+Mc 


+ 
RETURN 


( ): SUBROUTINE 


FIGURE D.3 OUTLINE OF PROGRAM FOR MOMENT-CURVATURE RELATIONSHIPS 
(NAPB) IN (MC) 
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(READ A) 
I5 (MAIN) 
1615 
16F5.0 
(READ B) 
1615 
16F5.0 


( ): SUBROUTINE 


FIGURE D.4 INPUT DATA 
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List of Computer Programs 


IMPLICIT REAL*8 (A-H,O-Z) 
COMNON/TEA2/SIGY,YS,AH,RSN,RES,FYP,EPM, SIGU 
COMNON/BL1/N,N1,AL, NM, NR(10) 
COMMON/BL2/CIP,CIN,AKP,AKN,ATP,ATN, BYP, BYN 
COMMON/BL3/AS,S1,X,CI,AK, BY, EQI, AT 
COMMON /BL5/SP (5) ,SCR(6,50,4) , PRO (4,100) 
COMNON/COS1/CU (6,61) 
COMMON /BL4/XX (4,100) , YY (3,100) 
DIMENSION BM (100),R(100),SR(100) , DEF (100) ,SD(100) 
DIMENSION 0Q(100),SH(100) ,P (100) ,BL(10) 
ELN=3.0*10.0*%4 
CALL READA 
DO 16 t=1, 100 
poe d=12 a 
sb 4d baler tig ELK) 
DO 12 d=1,2 
5 bl Cvinega  Retae) 
CONTINUE 
WRITE (6,5) 
WRITE (6, 4) 
CALL DEBU (1) 
CALL MP(1,BMU) 
CIP=EQI*ELN 
AKP=AK 
ATP=AT 
BYP=BY 
WRITE (6,5) 
WRITE (6, 3) 
CALL DEBU(~-1) 
CALL MP(-1,BMU) 
CIN=CI*ELM 
AKN=AK 
ATN=AT 
BYN=BY 
CALL MC 
READ (5,2) NL 
NNL=0 
CALL READB (P,Q) 
CALL DERM(Q,P,BM,SH,R, DEF, SR, SD) 
CALL CPRINT(SH,BM,R,DEF, SR, SD) 
Dot 150 
HO 6 =1,3 
JJ=J+3 
IF (J. EO. 1-.OR.J. EQ. 2) CU (JJ,1) =-SCR(6,1,0) 
TE (EOS) CU (Jd 1) 9Ch (6,.L,0) 
CONTINUE 
CUI=CU (4, 1) *SIGY/ (CU (5, 1) *CU (6, 1) ¥ELM) 
IF(NM.EQ.0)GO TO 14 
CALL NLE(BM, BL) 
DO 13 I=1,NM 
BS=BM (NR (I) ) 
BSS=DABS (BS) 
IF (BSS.LT. (CU1*1.05))GO TO 13 
PS=BSS/BL (I) 
CALL LOCAL (20, BL (1) ,CU1, PS, E1G1) 
CALL MX(BM,PRO) 
CONTINUE 
NNL=NNL+1 
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LE(RNL. LT. NL) GO TO 1 
CORLORS 


FORMAT (1015) 

FORMAT (/5X,'$$$$ NEGATIVE MOMENT REGION $$$$') 
FORMAT (/5X,'$$$$ POSITIVE MOMENT REGION £$$$$") 
FORMAT (/20X, '$$SE$$SFS$$S$SSH$$") 

FORMAT (/5X,8E14.5) 

STOP 

END 


SUBROUTINE CPRINT(SH,BM,R,DEF,SR,SD) 

IMPLICIT REAL*8 (A-H,O-Z) 

COMMON/BL1/N,N1,AL,NM,NR (10) 

DIMENSION BM(100),R(100),SR(100) , DEF (100) ,SD(100) ,SH(100) 
WRITE (6, 2) 

DORIS, N14 

WET TH (6993) 1 ,SH (1) ,BM (I) - B (1) , DEF (1) , SR (1) , SD(1) 

FORMAT (//5X,'**** BENDING MOMENT AND DEFORMATION ****', 

i Ape SHEAR, VOX, 9 EN", 18x, REG 1SY, UDBFUPISK ptSR! ,13%,15D") 
FORMAT (5X,13,2X,6E15.6) 

RETURN 

END 


SUBROUTINE MX (BM, PRO) 
IMPLICIT REAL*8 (A-H,0-Z) 
COMMON/BL1/N,N1,AL,NM,NR (10) 
COMMON/TEAY/BC, TC, FC, KT 
DIMENSION BM(100),PRO (4,100) 
NMI=NN+1 
WRITE (6,11) 
DO 1 I=1,NM1 
BE=0.0 
J2=1 
J3=N1 
IF (I. NE. 1) J2=NR (I-1) 
IF(1I.NE.NM1) J3=NR (I) 
DOMei=d2) 08 
IF (BM (J) «GT. BE) MB=J 
IF (BM (J) .GT. BE) BE=BM (J) 
CONTINUE 
SNM=PRO (2, MB) * (TC+ PRO (3, MB) -PRO (4, MB) ) 
WRITE (6,10) 1,SNM,MB 
CONTINUE 
RETURN 
FORMAT (/5X,'***##* MAXIMUM STRAIN ON CONCRETE SLAB ***#*%') 
FORNAT (/5X, '**#k',12,'TH SPAN IS',F7.5,'AT THE LOCATION’, I3) 
END 
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SUBROUTINE NLE (BM, BL) 
INPLICITREAL*8. (A-H, O~Z) 
CONMON/BL1/N,N1,AL,NM,NR (10) 
DIMENSION BL (10) ,BM(100) 
H=AL/DFELOAT (N) 

DO 1 I=1,NMN 

J=0 

NK=NR (I) -J-1 
TPN LE. 0).GO" "Tors 

IF (BM (NK)) 2,2,3 

J=J+1 

BL (LI) =DFLOAT (J) *H 

GO TO 4 

JJ=0 

NL=NR (I) +dd+1 

LF (NL.GE.N1)GO TO 6 
IF(BM (NL)) 5,5,6 

JIJ=JI5+1 

GO TO 14 

IF (JJ.GT.J) BL (1) =DFLOAT (JJ) *H 
CONTINUE 

RETURN 

HORM (/5X, 60 3) 

END 


SUBROUTINE READB(P,Q) 
IMPLICIT REAL*8 (A-H,0-Z) 
COMMON/BL1/N,N1,AL,NM,NR(10) 
DIMENSION P(100),K(16) ,A(16),B(10) ,Q(100) 
NM1=NM+1 

READ (5,1) (K (I) ,I1=1, 16) 
READ(5,2) (A(I) ,1=1, 16) 

READ (5,2) (B (1) ,1=1,-10) 

WRITE (6,5) 

DO 3 1=1,16 

UF (E CIys-20}080).G0 To 3 

WRITE (6,4) K (1) »A (I) 

CONTINUE 

WRITE (6,11) 

pO 10 §T=1,NN1 

WRITE (6,12) I,B (I) 
DOPMIE=A1,N1 


JIIZNA 
IF(1I.NE. 1) JJ=NR (1-1) 
IF (1. NE.NM1) JJJ=NR (I) 
PO (8 °IeJIfJdd 

Q (J) =B (I) 
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CONTINUE 

DO 9 I=1,NM 

QO (NR (I)) =O (NR (I)) 72.0 

Q(1)=0Q(1) 72.0 

Q(N1) =Q(N1) /2.0 

RETURN 

FORMAT (1615) 

FORMAT (16F5.0) 

FORMAT (1H1,5X,'**k*k LOAD POINTS ****'!_/) 
FORMAT (SX,*LOAD POINT AT',13,3X,F6.2,'KIPS."') 
FORMAT (/5X,'**** DISTRIBUTED LOADS ****") 


FORMAT (/5X,'****",73,'TH SPAN',F8.2,"KIPS/INCH ****!) 
END 


SUBROUTINE MP(K,BMU) 

IMPLICIT REAL*8 (A-H,0-Z) 
COMMON/TEA1/B,T,D,W,K1,K2,KF 
COMMON/TEA2/SIGY,YS,H,RSN,RES, FYP, EPM, SIGU 
COMMON/TEA3/RA,YRA,HR,RSNR,FYPR,DRS,RK 
COMMON/TEAY/BC,TC,FC,KT 

DW=D*W 

DM=DW*1. 1 

BT=B¥T 

BH1=D/2.0+T 

BA=DW+2.0*BT 

BTC=BC¥TC 

eel ba) 2 pee 

RAA=RA*YRA/SIGY 

Y4=D/2.0+RAA/W/2.0 

IF. (YU.. GEeD).GO TO 4 

ZB=B¥T*(T+D) +W* (Y4**2/2.0+ (D-Y4) **¥2/2.0) *1.1+ (D-Y4+T+DRS) *RAA 
BMU=ZB*SIGY 

WRITE (6, 10) BMU 

RETURN 

X=(BA-KAA) /2.0/B 

XX= =x 

Z=(XX+D+T/2.0) *BT+ (XX+D/2.0) *DM+ (XX*¥*2+X**2) *B/2.0+RAA* (DRS+X) 
BMU=Z*SIGY 

WRITE (6,10) BMU 

RETURN 

A3=B¥T* (D+T) /2. 

AU=W* DE*2Z/4, 

BMU=(A3*2.+A4) *SIGY 

WRITE (6,10) BMU 

RETURN 

Y¥2=+BA*SIGY/ (BC*0. 85*FC) 

BMU=BA*SIGY* (BH1+TC-Y2/2. 0) 

TRY cal ete aid TE (6, 10) RMU 

LF(Y2.LT.TC) RETURN 
RRS=(BA*SIGY-0.85*FC*BTIC) /2.0 

IF (RRS.GT. B*T) GO TO 20 

¥ 2=RRS/B 

BNU=SIGY* (B¥L* (D#tl+T/2.0-Y2) +D*¥W*1.1* (D/2.0+T-Y2) + 
1B* (Y2**2/2.0+ (T-Y¥2) **272. 0) ) +0.85*FC*BIC* (TC/2.0+Y¥2) 
WRITE (6,10) BMU 
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RETURN 
Y 2=(RRS-B*T) /W 
BMU=SIGY* (B¥T* (D+T) +W¥*¥1. 1* (Y2**2/2.0+4 (D-Y2) **2/2.0) ) 


1+0. 85*FC#BTC# (TC/2. 0+ T+Y2) 


WRITE (6,10) BMU 

RETURN 

FORMAT (10X,E15.7) 

FORMAT (5X,'**** IDEALIZED PLASTIC MOMENT IS',F9.1,' KIP-IN****') 
END 


SUBROUTINE DS(SH,PRO,SD) 
IMPLICIT REAL*8 (A-H,0-Z) 
COMMON/BL1/N,N1,AL,NM,NR (10) 
COMMON/BL2/CIP,CIN, AKP, AKN, ATP, ATN, BYP, BYN 
DIMENSION BM (100) ,SH(100) ,SD(100) ,S (10) ,SG (100) ,EM(100) 
DIMENSION PRO(4, 100) 
NM1=NM+1 
AH=AL/DFLOAT (N) 
DO 1.I=2,N 
FY=PRO (2,1) 
YP=PRO (3,1) 
Y D=PRO (4,1) 
Tevet. 060) Ke—1 
TE EY Cr OL 0) haul 
IF (FY. EQ.0.0) K=0 
TF(FY.LT. 0.0) CI=CIN 
TP(FicGh, 0. 0) CleclP 
IF(SH(I) .NE.0.0) BS=DABS((FY*CI-PRO(1,I))/SH(I) ) 
IF(SH(1I) .NE.0.0) FYM=FY*AH/BS 
IF (SH (I) .EQ.0.0) FYM=0.0 
CALI, SHERC(K, 0+ Us 0u.0,0.0,F YM, YP, YD, AH, SG (1)) 
PRCGHI) ON Ee Os 0) 56 (1) =s6 (1) oH (1) 
Die DFLOAT (1) 
CONTINUE 
SC ye oo (2) 
SG (N1) =SG (N) 
po. 18.%=1,N1 
SPi =U 
bos io el] 1,4 
alee 
AR uj 
IF (I. NE. 1) JJ=NR (I-1) 41 
IF (I.NE.NM1) JJJ=NR (I) -1 
Dom ede, od 
A=DFLOAT(JJJ-J) /DFLOAT (JJJ-JJ) 
MOm20eh= Jd, 0 
SD (J) =SD(J) +SG (K) *A 
Doe einkso dud 
SD (J) =SD (J) +SG (K) * (A-1.0) 
CONTINUE 
RETURN 
FORMAT (5X, 515) 
FORMAT (5X,8F15.5) 
END 
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SUBROUTINE LOCAL (N,AL,BMC,PY,EIG1) 
IMPLICIT REAL*8 (A-H,0-Z) 

COMMON/TEA1/B,T,D,W,K1,K2,KF 

COMMON/TEA2/SIGY,YS,AH,RSN,RES, FYP,EPM,SIGU 

DIMENSION H1(20) ,H2(20) ,H3(20) ,H4 (20) ,H5 (20) ,WA(20,20) ,WAA (20, 20) 

ELM=3.0*10. **4 

EIW=ELN*B¥*3*T**3/144 .0 

AFAC= (B¥*2+T**2) /12.0 

G=FLN/2.6 

GJID=G *B*T**3/3.0 

NO=N-1 

NT=N-2 

N1=N+1 

aca 

AMC=PY*AL- 

BL=AL* (1.0-BMC/AMC) 

IF (BL.LE.0.0) RETURN 

H=BL/DFLOAT (N 1) 

H10=H** (-2) 

H20=H 10**2 

DO 1=15N 

DIM=DFLOAT (I) /DFLOAT (N) 

H1 (I) =BMC+ (AMC-BMC) *DIM 

HM=H1 (I) 

CALL YAMA(HM,AF,Y1) 

O1=-AF*Y1 

O2=STRES (01,SIGY,FYP,AH, SIGU) 

CALL PORT (01,ER,GR) 

H2(1I) =EIW*ER *H20 

H3 (1) =GJD*GR¥*H10 

H4 (I) =B¥T*O 2*AFAC*H10 

H5 (I) =ELM*W**373.0/Y1 

CONTINUE 

DOs 2atl aay N 

DO ge el dan 

WA(1I,J) =0.0 

WAA (I,J) =0.0 

CONTINUE 

DO, 3, I=1,. 

I1=I-1 

LOANED 

TF (is EQ. 1).GO TO..4 

VEC lee Gls Nuc LO ye) 

Z1=H2 (11) +4. 0*H2 (I) +H2 (12) 

Z2=H3 (11) /2.0+H3 (I) +H3 (12) 72.0 

Z3=- (H4(I1) /2.0+HN4 (1) +H4 (12) 72.0) 

GO MLOR6 

Z 1=6.0*H2(1) +H2 (2) +0.5*H3 (1) 

Z2=1.5¥*H3 (1) +H3 (2) /2.0 

7 3=- (H4(1) +H4(2)) 72.0 

GO TO 6 

Z 1=+H2(NO) +6. 0*H2(N) +0.5*H3 (N) 

Z2=H3 (NO) /2.0+1.5*H3 (N) 

Z3=-2.0*H4 (N) 

{A (I,1I) =Z1+Z2+H5 (I) *.666667 
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WAAL, 1) = 22 

CONTINUE 

DORs (mele) iN 

2a 

bee Nee) tie (2) ) (os (1) +83 (12) ) /2.0 
Z2= (H4 (TL) +H4 (£2)) 72.0 

WA (12,1) =Z21+H5 (1) *. 1666667 
WA(1, 12) =WA (12,1) 
WAA(I2,1)=22 

WAAWIy 12) =22 

CONTINUE 

DiGmomeG= le Nl 

r2=14 1 

I3=I+2 

WeGleveel atl pa dae) 

WA (1,13) =WA (13,1) 

CONTINUE 

CALL SHUKI (WAA,WA,N, EIG1, 4) 
RETURN 

BORMAT (/5%7 5.5) 
BORMAN Zo Xap Ee) 

END 


SUBROUTINE YAMA (AM,AF,Y1) 
IMPLICIT REAL*8 (A-H, O-Z) 
COMMON/COS1/CU (6,61) 

MT=-1 

IF(HT.EQ.1) MRI=1 

IF (MT. FQ. 1) M0=2 

TF (MT.EQ.1) MK=3 

IF (UT. EFQ.-1) MRI=4 
IF(MZ.EQ.-1) MO=5 

IF (MT. FQ.-1) MK=6 

HAM=CU (MRI, 1) 

IF (DABS (AM) . LE. DABS (HAM))GO TO 1 
MW=1 

MW=MW+1 

H1=CU (MRI,MW) 

IF (DABS (AM) «GT. DABS (H1)) GO TO 2 
MW1=NU~1 

H2=AM-CU (MRI, MW1) 

H3=H1-CU (MRI,MW1) 
H4=CU(MO,MW) -CU (MO, MW1) 
AF=CU(MO,MW) +H2*H4/H3 

Y¥1=CU (UK, MW) +H2* (CU (MK, MW) -CU (MK,MW1) ) /H3 
RETUKN 

AF=CU (MO, 1) *AM/HAM 

Y 1=CU (NK, 1) 

RETURN 

END 
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SUBROUTINE PORT (EP, ER, GR) 
IMPLICIT REAL*¥8 (A-H,O-Z) 
COMMON/TEA2/SIGY,YS,H,RSN, RES, FYP, EPM, SIGU 
ELM=3.0*10.0**4 
EP=DABS (EP) 

A=EP-FYP 
THER) AW, SIGIZ 

1 BRE wD 
GRart.0 
RETURN 

2 EST=ELN/H 
EPST=RSN*SIGY/ELM 
AA=FST* (0. 2-EPST) 
A1=S1GU**2-SIGY** 2-2. 0O*¥SIGY*AR 
A2=2.0* (SIGU-SIGY-AA) 

A3=A1/A2 

B1= (SIGY-A3) * (SIGY-A3-2.0*EST*EPST) 
C1=2. 0 *EST* (SIGY-A3) 

ETAN=C1/DSQRT (B1+C1*EP) /2.0 
ER=BTAN/ELM 

GRE 2.16 £01,344. 07 ER 4.0) 

RETURN 

END 


SUBROUTINE SHUKI (DD,CC,N,RMN1) 
Dee eR EAL OWA Hy OZ) 
Bigpis Lom Z.(20)",¥ (20), YY (20) - H (20720), C (20,7 20) 7D(20, 20), HH (20,20) 
1,CC (N,N) ,DD(N,N) 

DO 120 I=1,N 
DO 110 K=1,N 
Gu7h) =CC (1, K) 
Dak). =DD (17K) 

110 CONTINUE 

120 CONTINUE 
He(HRTEO. 1) GO TO 210 
GHG LORY OM NTC) 
DO 140 I=1,N 
DO 140 K=1,N 
ABC=0.0 
DO 150 J=1,N 
ABC=ABC+tC (I,J) *D (J, K) 

150 CONTINUE 
HB (I,K)=ABC 
H(1I,K) =ABC 

140 CONTINUE 
KK=0 
NA=N 

10 2Oe 207 t= 1 NA 
¥ (I) =1.0 

20 CONTINUE 
KK=KK+1 
KA=2* (KK/2) 
HID (AE SG TNO) S G/N), Y (NA) =-1.0 
B=0.000001 
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40 DO 50 J=1,NA 
YY (T) =0.0 
DO 60 J=1,NA 
YY(T) =¥Y (1) +H'(1,0) *¥ (0) 

60 CONTINUE 
IP (ree al) A=YY (1) 

_ 50 CONTINUE 

PECUADS (A) * Gl. 10, UE Io) 8GO TO 56 
DO 55 I=1,NA 
Y¥ (I) =Y (1) +DFLOAT (1) 

55 CONTINUE 
GO TO 40 

56 DO 70 I=1,NA 
NE ((5E)) SNES (GE) PATA 

70 CONTINUE 
S=A/B 
B=A 
IS» ols Ws OOOO) Eo Wo Le 
WES» ollie Wa DII9S)). EOr ze We 
1 (OSS) BlHOs- Na) Or GO. 7S 


AA=A 
DO 71 I=1,NA 
ZT) =v (1) 
71 CONTINUE 
GO TONTO 
75 TP(Ae: GEL AR) GO TO’ 80 
A=AA 
DO 76 I=1,NA 
YU) 221) 


76 CONTINUE 
SOV ITIRRe Fon 2) GO) TOs 200 
PAW ny Nasi 
60"'T0"220 
210 AED(45 13/01, 1) 
220 WRITE(6,230) A 
230 FORMAT (/1X,'THE FIRST EIGEN VALUE IS',E12.4) 
RM 1=A 
IF(N .EQ. 1) RETURN 
WRITF(6,240) (Y¥(1) ,1=1,N) 
240 FORMAT ( 3X,"CORRESPONDING MODE **'/(3X,8E14.4)) 
RETURN 
END 
SUBROUTINE TOKYO (NO,A) 
IMPLICIT REAL*8 (A-H,0-Z) 
DINENSION AC(20) ,AB(20) ,A (20,20) 
NO1=NO-1 
Rilo 1) = 107k (let) 
IF(NO .EQ. 1) RETURN 
DO 80 N=1, NO1 
DO 50 I=1,N 
AB(I) =0.0 
Ret 0.0 
pov50 J=1,N 
AB(I) =AB(I) +A (I,d) *A (3, N+1) 
AC (I) =AC(I) +A (N+1,J) *A (J, 1) 
50 CONTINUE 
ACB=0.0 
DO 60 I=1,N 
ACB=ACB+AC (I) *A (I, N+1) 
60 CONTINUE 
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A (N+1,N+1)=1.0 +4),N4+1)-ACB 
DO 70 TN a eee eas © 
A(N+1,1) =-A (N+1,N41) ¥AC (T) 

A (I,N4+1)=-AB (I) *A (N41, N41) 
CONTINUE 

DO 80 I=1,N 

DO 80 J=1,N 

A (I,J) =A (I,J) -A (I,N+1) *AC (J) 
CONTINUE 

RETURN 

END 


SUBROUTINE DEBU (K) 

IMPLICIT REAL*8 (A-H,O-Z) 
COMMON/BL1/N,N1,AL,NM,NR (10) 

COMMON PES/AS ,S1,X,Cl, Ak, BY, EO1,LA 
COMMON/TEA1/B,T,DD,W,K1,K2,KF 
COMMON/TEA2/SIGY,YS,H,RSN,RES, FYP, EPM,SIGU 
COMMON/TEA3/AR, YRA,HR,RSNR,FYPR,DS,RK 
COMMON/TEAY/BC,TC,FC,EN,KT 

FORMAT ( /5X,!**** COMPUTATIONAL RESULTS ****', /5X, 
1 YaREA COR STEED SBEAM (EN) ',FS8.3,5X, 

2"MOMENT OF INERTA OF STEEL BEAM (IN**4)',F8.3,/5X, 
3*NUETRAT AXIS “OF COMPOSITE BEAM (IN)',F8.3,5X, 
& "MOMENT OF INERTIA OF COMPOSITE BEAM',F9.3,/5X, 
5'*SHEAR CORECTION FACTOR OF COMPOSITE BEAM',F8.3 ) 
D=DD+2.0*T 

SLL=AL/2.0/DFLOAT (NM+1) 

Ev(K. EQ. 1) GO TO 7 

IF (K.EQ.-1)GO TO 8 

A1=BC/EN 

RZ2=TC 

GO TO 9 

Al=DS 

A2=AR*RK 

AS=2.0*B*T+ (D-2.0*T) *W 

GATT, INU (Ke NS AD) iN fA.2)7-X) 

GALL MISE (BT PD, W, ST) 

GPU MISG(K a, Dp oie, AS, Al, AZ, TA,CL) 

CAIN ISCR (K,5,T,;D,W,Xehl,A2; Ci; AK) 

Whim (6, 10) AS, SI,X,C1 AK 

De (Hob O..1). CALS BOQUD(AS, SLL, EN; BC, TC, D;SIGY,SL,4, EQL) 
IF.(K. EQ. 1) SY=SIGY*EQI/ (D-2Z) 

IF (K.EQ.-1) BY=SIGY*CI/ (D-X) 

WRITE (6,11) BY 

FORMAT (5X, "YIELD MOMENT',F8.1) 

RETURN 

END 


SUBROUTINE NUE (K,AS,D,A1,A2,X) 


rds | 
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(1.8%, "ReHOM ase 


ae he 


IMPLICIT REAL*8 (A-H,O-Z) 
TE(ReECc dy co LO 4 

IF (K.EQ.-1)GO TO 2 

BN=A1 

Tea i? 

AN=AS*D/2. 0-BN*¥TC¥*2/2 .0 
IF(AN.LT.0.0)GO TO 3 
X=AN/ (BN*TC+AS) 

RETURN 

AA=2.0*AS/BN 

BB=D/2.04TC 
CC=AA**24+4.0*AA*BB 

X= (-AA+DSQRT (CC))/2.0-TC 
RETURN 

AR=A2 

DS=A1 

X= (AS*D/2.0-AR*DS) / (AS+AR) 
RETURN 

END 


SUBROUTINE MISB(B,T,D,W,S1) 
LHPLICLT REAL*S@ (A—H, 0-7.) 
DD=D-2.0*f 

A1=B¥T**3/12.0 
A2=W*¥DD**3/12.0 

A3Z=B¥T* (DD/2.0+4+T/2.0) **2 
ST=2.0*A1+A2+A3*2.0 

RETUPN 

END 


SUBROUTINE MISC (K,X,D,SI,AS,A1,A2,TA,CI) 
IMPLICIT REAL*8 (A-H,O-Z) 

Peirce 1) GONTOM 1 

IF(K.EQ.-1)GO TO 2 

BN=A1 

TC=h2 

Tr (2 LT.0.0)GO.TO 3 

TA=BN*TC+AS 
AI=BN*TC¥* 3/12. 0+ BN¥TC¥ (X+A2-TC/2.0) **2 
CI=SI+AI+AS* (D/2.0-X) *¥*2 

RETURN 

TC=TC+X 

GO TO 4 

DS=A1 

AR=A2 

TA=AS+AR 
CI=SI+AR* (X+DS) **2+AS* (D/2.0-X) **2 
RETURN 

END 
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SUBROUTINE SCF(K,B,T,D,W,X,A1,A2,AI,AK) 

IMPLICIT REAL*8 (A-H,0-Z) 

DD=D-2.0*T 

AS=2.0*B¥T+D*U 

TECK. FO. 1) Go T0604 

IF (K.EQ.-1)GO TO 2 

BN=A1 

TC=A2 

IF(X.EQ.0.0) TC=TC+X 

AT=AS+ BN*TC 

G1=-BN/2.0 

EX=TC+ X 

) BF 663.0000) eT C 
G3=-TX**2+61 

R= x ; 

GE (ROL TA0 10), TL=On 0 

GALL? OX2(G61,\00;, G3, TX, TT, 001) 

QQ=Q001/BN 

GO=BN * (TX**2-TT**2) 72.0 

G1=-B/2.0 

G3=GO+X**2*B/2.0 

X1=X-T 

GalAION2 (G1, 0.0,63, X, X1,001) 

Q0=00+001/8 

GO=GQtB* (X**2-X1**2) /2.0 

G1=-W/2.0 

G3=GQ+X1**2*4/2.0 

X2=X1-DD 

GhULTION2 (G 90,10 (68 HC, X2 007) 

QQ=00+001/¥ ; 

GQ=GQtW* (K1*%*2-X2**2) /2.0 

Gi=-B/2.0 

G3=GQ+X2** 2*B/2. 0 

X3=X2-1 

GhtIgION2 (61,02 0,63, 12, %2, 007) 

Q0=00+001/B 

Gouro; 5 

DS=A1 

AR=A2 

AT=AS+AR 

GQ=AR* (X+DS) 

Q0=0.0 

GiOy TO. 

AK=AT*QQ/AI¥¥*2 

RETURN 

END 


SUBROUTINE QX2(A,B,C,D1,D2,ATL) 
IMPLICIT REAL*8 (A-H,0-Z) 


INTEGRATION OF (AX**2+BX+C) **2 BETWEEN D1 AND D2 
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RIS RAE 2 4 (D1 eS D2 4*S)07 5 20 

A2=A*B* (D1*¥*4-D2¥*4) 72.0 

A3= (B¥*¥*2+2.0*A*C) * (D1**3-D2**3) 73.0 
A4Y=B*¥C* (D1**2-D2**2) 

AS=C*¥* 2* (D1-D2) 

AIT=A1+A2+A3+A4+A5 

RETURN 

END 


SUBROUTINE EQUI (AS,AL, EN, BC,TC, D,SIGY,AL,DX,EQ1) 
IMPLICIT REAL*8 (A-H,O-Z) 

BN=BC/EN 

ELM=3.0*10.0**4 

F1=AS*SIGY 

F2=4.0*BC#TC 

IF (F1.GT.F 2) F=P2 

IF (F1. LE. F2) F=F1 

AK=ELM*BN/63.6/F/AL 

AP=EN* AS/BC 
SS=1.04+(D+2.0*TC) *(1. O+AK*AP) /AP 

Z=AP* (DSQRT (SS) -1.0)/(1.0+AP*AK) 
TE(ZIGT.1C) GO 10 42 

DY=AK¥*Z**2 

EQI=AI+AS* (D/2.0+TC-Z-DY) **2 

EQ 2=BN *Z** 3/12. 0+ BN¥Z**2¥* (DY+Z/2.0) /2.0 
EQI=EQ1+EQ2 

K=TC-7 

DX=DY-Z 

WRITE (6,1) AK,X,DY,EQI 

FORMAT (//5X,'**** EFECT OF SLIP ****'"',//5X, 


T*AK (=DY/2"% Dy" PFE. 3,5K;, 

2'NUETRAL AXIS FRCM TOP OF STEEL BEAN',F8.2,/5X 
3,"DISTANCE OF TWO NUETRAL AXIS',F8.2,5X, 

W* EQUIVALENT MOMENT OF INERTIA (IN¥#*4)',F12.2 ) 


RETURN 

AK 1=2. O*AK*TC 

AK2=1.0+AK1 

AK3=1. 0+AK1/2.0 

AC=TC*BN 

PZ=AS* (D/2.0+AK 3¥*97C) +AC*TC/2.0 
Z=PZ/ (AS*AK2+AC) 
DY=AK 1% (Z-TC/2.0) 
EQI=AI+AS* (D/2. O-AK2*Z+AK3¥*1C) ¥*2 
FQ2=BN*TC#*3/12.0+AC* (Z-TC/2. 0) *AK2* (Z-TC/2.0) 
EQI=EQ1#EQ2 

GO TO 3 

END 


SUBROUTINE RADA 
IMPLICIT REAL*8 (A~H,O- 2) 


DT#Da +0. 
STA (ST. TH. 
SN 


Casas) AAD, i aei te seat 


$*¥ (¥a+S-3T+0.¢ a) «; 
0.S\ (0. S\S+¥ ad) #S*HS THEO. SENE 


yd 5 
,XEN\\, ee ew BTTe IO TOM pa ae 
XO cE OE, * conte 
K2\,S- BF, "MART JIZST2 10 SOT HAY SIA 
.¥2,S.8%,"22NA JABTAUH ORT tO 4 
{ $.Sta,’ (MeeHT) ALDATHT 7O THIHON 


(0. iseceie ete <oP-k) #9490,5 


fae) 


COMMON/TEA1/B,T,DD,W,K1,K2,KF 
COMMON TEA 2Z/SLGY,YS,H,RSN,RES,;FYP,EPM,SIGU 
COMMON/STEA3/RA, YRA,HR,RSNR,FYPR,DRS,RK 
COMMON/TEA4Y/BC,TC,FC, EN, KT 
COMMON/BL1/N,N1,AL,NM,NR (10) 
READ(6,1)8,1,D,¥ 
FEAD (57) 00 Wal, RSN,SIGU 
READ (5,5) K1,K2,KF 
READ(5,1)RA,YRA,DRS,HR,RSNR, RK 
READE (oro) ps Gr LG yh Ge EN eKce 
READ OPM ls, NyiNM PONE (L) t= 11,10) 
DD=D-2.0*T 
N1=N+1 
ARIE (GO GZ) ln PS i i 
SEAS, (SUD) Sine 5 al HRS SIC 
WRITE (6,8) RA, YRA,DRS, HR, RSNR, RK 
WEEE (Gal 0) SGC iG iN 
_ WRITE (6,12) AL,N,NM, (NR(I) ,I=1, NM) 
1 FORMAT (8F6. 0) 
2 FORMAT (1H1,//5X,'**** STEEL BEAM DIMENSION ***x!"t, 
W//5x, FLANGE OWIDTH"' (F843 5X ,"FLANGE THICKNESS=',F8.3, 
2/5X, ‘BEAMHEEGHT=!,F8.3,5X,'WEB THICKNESS=# , F8. 3) 
H FORMAT (//5X,'**** STEEL PROPERTIES ****!, 
J foie Le Dont Gos Oar Oks RAL OOF Ly Biol) t,t Bs 2, 
2fOL RAL LOLOT REPS Ao Ly BES (Ht phoaey OX, DULLINAT EH eoOTR ESS =" 7 Fo. 3) 
8 FORMAT (//5X,'**** REINFORCEMENT ****?! , 
Woe UR A= op NOMA pou, Y Ei D SUR HS S— lp F Ole 2k, | D LowwANG =i Py Oriy 
2f OX pL LOOT E/E Sa) = EOw2ig kept Ra lLOwOFreEPoNS iL) / GES (Ye boc, 
Bia kage eels Galak N Gea te TO 3) 
10 FORMAT (//5X, '**** SLAB PROPERTY ****%'", 
UA D068 5 OSTWING WALDEE TIS” 516) 5 Wp Shp USIOING SCARS SaY pio5 Zp D5 
2/5X 7 CONGRETERSTRENGTH=" )FO.2, 5X, 'N=E(S) /E(C)=",F 6c) 
A VONRSUINIE (ZS 5 U SSEeE TTY T TEAMED S| facets U 
1/75xk, CTOPAE BEAN LENGTH=",F8.2,5X, "NUMBER OF SEGHENTS<* 13, 
2/5X, "NUMBER OF INTERMEDIATE SUPPORTS=!,13,' AND THE NUMBERS OF 27 
3'THE INTERMEDIATE SUPPORTS=!,1015) 
SeFORM AL (Ou) 
9 FORMAT (4F6.0,15) 
47 FORMAT (Fo. 041215) 
RETURN 
END 


SUBROUTINE REACT (M,0,P,EM,BM,EI,G,X) 
IMPLICIT REAL*8 (A-H,0-Z) 
COMMON/BL1/N ,N1,AL,NM,NR (10) 
DIMENSION P(100),EI(100) ,X(10) ,S(10, 10) ,Q(100) ,4 (4) ,G (100) 
DINENSION BN (100),SH(100) ,R(100) ,EM(100) ,DEF(100), 
1PP(100) ,M1(4) , TM (100) 
CALL FDM(M,Q,P,EM, BM, EL,G,SH,R, DEF) 
DO 1 J=1,NM 
1 ¥(J)=-DEF (NR (J) ) 
poe pte (.8 
4 M1(I)=0 
M1(2)=1 
DO 2 J=1,NM 


‘MH, Anes a<a 2 
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8-832," (2 peya\3 $0 OLTAA'.AC,5, 84 '=geuat2 ‘\ 
(£.8%, *=223HT2 @DANTT IU! X2,o.845 * EY) 294\, (C2) 284 FO OLTAS! .LANS 
, ree PERRIDANINISH see? X2\A)T 8 
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a, evoug.1o12a" . 
heey YOREIONS ange wees SENN) TAMHOT ‘or. 

2.8.94, '=BAMWADIND BhIe' Xe .f. @s,'=R?arw gAl2’ x 
(S.97,* - (oy sy (2) a=" teh. 31, =n vowRare ATSAIUOD * Xe 
.'eer* HTONSI NAT #494! X2\\) TANROT sr 
,€I\*=eTuahode To asanuh? xe (2.89, Hh TORI L BARA JATOT’ XO\\T 
.* TO 2f24H0K BAT GKA ',ET, "seTROUgveE SATCANAATHE FO SIAHUH" Xe 
(210t, '=2TsO09402 TLALGMNASTYT BRTIE 
{2T0r) TAMRO’ 2 
(21,0.999) TANSOT @ 

(2S ,0.07) octg rt 


(X,, 18 nei .M) T2AgR giro 
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DO 5 I=1,N1 

PP (I) =0.0 

IF (I.EQ.NR (J) ) PP (1) =1.0 
CONTINUE 

CALL FDM(M1,TM,PP,TM,BM,EL,G,SH,R,DEF) 
POalor tat at) 

BM (I) =BM (I) -EM (1) 

DO 3 K=1, NM 

S (J, K) =DEF (NR (K)) 
CONTINUE 

CALL SMA(S,X,10,NM) 
RETURN 

END 


SUB OUR INE aE DINGO Oly Ely full gS hip Ey le piGrgiarg cv y) DE) 
AAVME AES ITE WN Posts) (Use) 8 0) 7A 
COMMON/BLI/N,N1,AL, NM, NR(10) 
DIMENSION Q (100) ,P (100) ,8 (100,100) ,C (3) ,BM (100) , EI (100), 
TATA U EO BCE) Ope (00) Bee too) 
,EM (100) ,G (100) , PA (100) 
“hn =DFLOAT (N) 
=AL/AN 
thei 
DO) iS} Tle a 
Tr(M (1) - BO. 0) 0 (1) =0.0 
TLE (au (2M) 5 OG) )) 12 (46) S05 © 
I (HL CS)) os KOS O) ISI) Oe 
IF (M (4). EQ. 0) DEF (1) =0 me 
BM (I) =0.0 
CONTINUE 
HN=H/6.0 
DO 1 Jei,N1 
AI=DFLOAT(I-1) 
DO 1 J=1,N1 
AJ= le 1) 
JJ=I- 
eh Sete 
B (1,J) =H*AI* (AN-AJ) /AN 
GOnLOmat 
B (1,d) =H*AJ* (AN-AI) /AN 
CONTINUE 
C (1) =HN 
C (2) =4.0*HN 
C (3) =HN 
DO mes aei—aley Nat 
PA (J) =P (J) 
Ova /mea—alyrs 
Tpdss slit liee, 
Ltd .80.0.0R.id- EO. (N1*1)} GO TO 7 
BA CePA (i) +0 (0) CT) 
CONTINUE 
DO 9 J=1,N1 
po 9 J=1,N1 
Bhi (J) =BM (J) +PA (I) *B (I,J) 
DO m2 On =a) Niat 
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IF(IJ.EQ.0.OR. IJ. EQ. (N1+1))GO TO 23 
Q01(J) =001(d) +BM (Id) /EI (1d) *C (I) 
CONTINUE 

CONTINUE 

DO’. 34 'T=1,'N1 

DO 34 J=1,N1 

DEF (I) =DFF (I) +B (I,J) *Q01(J) 

DO". E41 N 

S1=BM (I) /EI (I) 

S2=BM (I+1) /EI (I+1) 

BE1=BM (I) -EM (1) 

BE2=BM (i+ 1) -EM (I+1) 

S (I) = (EE2-BE1) /H+ (2.0*Q(I) +Q(1+1) ) *HN 
R (I) = (DEF (I+1) -DEF (I) ) /H+ (2.0*S1+S2) *HN 
CONTINUE 

BE1=BM (N) -EM (N) 

BE2=BM (N1) -EM (N1) 

S (N1) = (BE2-BE1) /H- (2.0*0(N1)+0(N)) *HN 
R (N1) = (DEF (N1) -DEF (N) ) /H- (2. 0*2BM (N1) /EI (N1) +BM(N) /EI (N) ) ¥HN 
RETURN 

END 


SUBROUTINE NAPB(FY,BM,Y1,YD,SL) 
IMPLICIT REAL*8 (A-H,O-Z) 
COMMON/TEA1/B,T,D,W,K1,K2,KF 
COMMON/TEA2/SIGY,SY,H, RSN, RES, FYP, EPM, SIGU 
COMMON/TEA3/RA, YRA, HR, RSNR,FYPR,DRS,RK 
COMMON/TEAU/BC,TC,FC,EN,KT 
DIMENSION GB(100) ,cY(20),¥ (30) 
BMM=BM 

ELM=3.0*10.0**4 

FYP=RSN*SIGY/ELM 

FYPR=RSNR¥YRA/ELM 

AS=2.0¥*T*B+W*D 

IF(FY.LE.0.0)MT=-1 
Prthy Gt. 0.0) MT=1 

PYT=20. 0 

BUI=0.0 

MR=0 

T1=L/FLOAT (K1) 

D1=D/FLOAT (K2) 

T2=T1/2.0 

D2=D1/2.0 

IR=1 

K5= (K141)* (KF+1) 

KG= 2* (K2+1) 

AT=2.0*B*T+W*D 

Ce=000 

Y¥2=T/2.0-¥1 

CALL FORCE(FY,Y¥2,B,T,KF,K1,1,GB,BDM) 
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BM=BDM 

NO 37 °1=17 6.5 

CC=CC+GRB (TI) 

Y2=T+D/2.0-Y1 

CALI: PORCH (PY, Yate WetDetl hoy te Gh, DDN) 
BM=BM+BDM 

DO 32 I=1,K6 

CC=CC+GB (1) 

Y2=T+T/2.0+D-Y1 

CAUINSFORCE (FY, Y2,B,T, KF,K1,1,68,5BD4) 
BM=BM+BDM 

DOV 35 lS (KS 

CC=CC+GB (I) 

aCe (Uelae)) WE ashy as) 

SK=FY* (D+2.0*T-¥14DRS) *RK 
CPr=STRES(SK,YRA,FYPR,HR, SIGU) *RA 
CC=CCt+CF 
BU=BNtCHe (De2. OF TY 14+ DRS) 
AT=AT+RA 

GOVTO Ste! 

PT=AS*SIGY 

F2=0.85*FC* BC*¥TC 

cA Ge ci Nee! GaP 8 ia PA 

UU=1.25*F1 

GO TO 22 

UUSP 2 5 
¥S=SLIP (DABS (CC) ,UU) /SL 
Y2=-TC/2.0-Y1+YS/FY 

YD=YS/FY 

CAE ORGHICE Yi Ni2 fb Opie WprKiL pe, Gb EDIT) 
BM=BM+BDHM 

KT2=2%* (KT+1) 

DOS Ger ik 2 

CC=CC+GB (I) 

AT=AT+BC*BT/10.0 

Y (IR) =¥1 

CY (IR) =CC 

CAT=CC/ (AS *SIGY) 

TF (DABS (CAT) .LE. (1./(10.**4) ) ) RETURN 
DP (DR. NEL A) .COCTO"S 

IF (CC)6,6,7 

Y1=0.9'8*yY 1 

IR=1IR+1 

GO TO 8 

YU=4, O24 V4 

IR=IR+1 

GO TO 8 

JR=IR~1 

LF (IR. GT. 20) RETURN 

EVYX=CYV(IR) —CY(IR) 

TPUDABS (CYX) LE. 10. 0** (~3) ) RETURN 
CYY=(" (ih) ) (OR) ) 7 (Cr (OR) Cy (LR) 
Vaal = Vo (audit Ca van(clR)) eG Yay 

JR=IR+1 

GO” TOs 

END 
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SUBROUTINE FORCE(FY,Y1,A1,A2,L1,L2,MA,GB,BM) 
ae Tee Cee: REAL*8 (A-H, 0-2) 
COMMON/TEA2/SIGY,YS,H,RSN,RES,FYP,EPM,SIGU 
COMMON/TEA4Y/ BC,TC,FC,EN,KT 
DIMENSION G(10,10) ,GA (3) ,GB (100) 

100 FORMAT (/5X,8F 15.5) 

L11=1.1+1 
L21=L2+1 
LL=L11*L21 
A11=A1/DFLOAT (L1) 
A21=A2/DFLOAT (L2) 
BM=0.0 
WO) 4! EES ey 
DO 1 J=1,1L11 
G(I,J3) =FY* (Y1+A2/2.0-A21*DFLOAT (I-1) ) 
G1=G (I,J) 
IF (MA.EQ.1)G (I,J) =STRES(G1,SIGY,FYP,H,SIGU) 
IF (MA. EQ. 2) G (I,J) =FOC (G1, FC) 
1 CONTINUE 
DO) Sy WS yb ind 2 
TRE" (ie E Orel ed = 
TENTSE OC L2H) TAS 
DOMGeOR— 1, UL Ile 
Toad. BO. .d) oc i=2 
Tr (@. £0<L Hj VAS LA 
ITJ=J+ (I-1) *L11 
GB(IJ) =4.0*G (I,J) +2.0* (G (11,3) +G (I,J1))+G(1i,J1) 
BM=BMN+GB (IJ) *(Y1+A2/2.0-A21*DFLOAT (I-1) ) 
YO=Y1+A2/2.0-A21*DFLOAT (I-1) 
6 CONTINUE 
PERO. 1) GOP TORS 
DOM S=H=2, ih 1 
Jt=J-1 
J2=J3+1 
Id=dt(2= 1) 4114 
GB(IJd) =8.0*G (I,J) +2.0* (G(I,J1) +G (I,J2)+G(I1,9))+G(I1,32) +G (I1,J1) 
1+2.0*G(1I1,J) 
BM=BM+GB (IJ) * (Y1+A2/2.0-A21*DFLOAT (I-1) ) 
YO=Y1+A2/2.0-A21*DFLOAT (I-1) 
7 CONTINUE 
5 CONTINUE 
DEFUL2: EQ 4) GO TOLVZ 
DO 8 I=2,L2 
r=1—1 
I2=1+1 
ne) 8) wes iby is 1 
IF(J.EQ.L11) J1=L1 
TPWOVEOLA) S1=2 
IJ=J+ (1-1) *L11 
GB (IJ) =8.0*G (I,J) +2.0* (G(I,J1) +G(I1,J) +G(12,J))+G(11,J1) +G (12,J1) 
1+2.0*G(I,d1) 
BM=BM+GB(IJd) *(Y1+A2/2. 0-A21*DFLOAT (I-1) ) 
YO=Y1+A2/2.0-A21*DFLOAT (I-1) 
9 CONTINUE 
LFi(L1..E0. 1) GO\TOL 8 
po 10 J=2,1L1 
. Id=J+ (I- 1) *¥L11 
Jt=J—1 
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V2 
cal 


J2=J3+1 


GB (IJ) =16.0*G (1,4) +4.0* (G (1,31) +6 (I,J 2) +G (11,5) +6 (12,J)) 


146 (11,91) +4G(12, 31) +6 (11,02) +G (12,32) 
BN=BN+GB(IJ) *(¥1+A2/2. 0-A21*DFLOAT (I-1) ) 
YO=¥14+A2/2.0-A21*DFLOAT (I-1) 

CONTINUE 

CONTINUE 

Dx Meoa feet a 

GB (1) =GB (I) *A11*A21/36.0 
BM=BM*A11*A21/36.0 
RETURN 

END 


FUNCTION FOC (SN,FC) 

IMPLICIT -REAL*8 (A-H,O-Z) 
Tio, Gre 0) FOC=0.0 

EP=-2.5%10.** (-3) 

E=2.0*DABS (FC/EP) 

IF (SN.LE.0.0) FOC=E*SN/ (1. 0+ (SN/EP) ** 2) 
RETURN 

END 

FUNCTION STRES(EP,SIGY,EPST,H,SIGU) 
IMPLICIT REAL*8 (A-H,O-Z) 
ELM=3.0*10.0**4 

EPY=SIGY/ELM 

EST=ELM/H 

EAI=EST**2/ (SIGU-SIGY) /4.0 
EA2=EPST+2.0* (SIGU-SIGY) /EST 
SIG=-EA1* (DABS (EP) ~EA2) **2+SIGU 

IF (DABS(EP) .LE.EPY) STRES=ELM*EP 

IF (EP.GE.EPY.AND.EP.LE.EPST) STRES=SIGY 
IF (DABS (EP) . GE. EPST) STRES=SIG*EP/DABS (EP) 
TF(EP.LE.-EPY.AND.EP.GE.-EPST) STRES=-SIGY 
RETURN 

END 

FUNCYION SLIP (P,USC) 

IMPLICIT REAL*8 (A-H,O-2Z) 

P1=DABS (P) /USC 

LF(P1.LT. 1.0) SLIP=-0.015189/ (P1-1. 11873) -0.013577 
TF (P1.GT..1..0).SLIP=P1-0.9 

RETURN 

END 

FUNCTION GAM (AA) 

IMPLICIT REAL*8 (A-H,0-Z) 
COMMON/TEA2/SIGY,YS,H,RSN,RES,FYP,EPN,SIGU 
ELM=3.0*10.0**4 

A=DABS (AA) 

G=ELN/2.6 

GAN=AA/G 

Se leat 2 

Gi =i Zc 

IF(A.GT. TY) GAM= ((A-TY) *H/G/4.0+GY) *A/AA 
RETURN 

END 
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SUBROUTINE SSP (SL,BM) 
IMPLICIT: REAL*8 (A-H, 0-2) 
COMMON/BLI/N,N1,AL,NM,NR(10) 
DIMENSION BM(100),S1L(10) 
NM1=NMN+1 

DO 1 T=1,NM1 
IF(I.NE. 1) JJ=NR (I-1) 

4039 (GES OG Dy yas 1 
IF(I.NE.NM 1) JOJ=NR (T) 

IF (I.EQ.NM1) JJJ=N1 

K=0 

DO 2 J=JJ,JJJd 

IF (BM (J) .GE.0.0) K=K+1 

SL(I) =DFLOAT (K) *AL/DFLOAT (N) /2.0 
CONTINUE 

RETURN 

END 


SUBROUTINE SHERG(K,FY,Y,YD,FYF,YF,YDF,AH,SG) 
IMPLICIT REAL*8 (A~H,O-Z) 
COMMON/TEA1/B,T,DD,W,K1,K2,KF 
COMMON/TEA2/SIGY,YS,H,RSN,RES,FYP,EPM,SIGU 
CONMON/TEA3/R2A,YRA, HR, RSNR, FYPR, DRS, RK 
COMMON/TEA4/BC,TC,FC,EN,KT 

DIMENSION GB(100) ,GC(100) ,S(11) 
TP(KjMNg2),3 

SK=-FY* (DRS+Y) *RK 

SKF=-FYF* (DRS+YP) *RK 

C=STRES (SK, YRA,FYPR,HR, SIGU) ~STRES (SKF, YRA,FYPR,HR,SIGU) | 
C=C#RA 

COMO” h 

C=ON0 

GO TO 4 

YA=-TC/2.0+YD-Y 

Y B=-TC/2.0+YDF-YF 

CHEL FORCE (FY, YA, BC, TC, 1,KT,2,GB, BDN) 
CALL FORCE(FYF,YB,BC,TC,1,KT,2,GC,BDM) 
KET=2* (KT+1) 

C2080 

DOl Ghia}, RTT 

C=C+GB (I) -GC (I) 

YA=Y-T/2.0 

YR=YF-T/2.0 

CALL FORCE (-FY,YA,B,T,KF,K1,1,GB, BDM) 
CALI, FORCE (-FYF,YB,B,T,KF,K1,1,GC,BDM) 
KTT= (K14#1) * (KF+1) 

Pom? I=1,KTT 

C=C+GB (I) -GC (I) 

DH=DD/10.0 

DA=DH*W 
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S (1) =C/W/AH 
DOSS eS 214 
YA=Y-T-DH*D FLOAT (1-1) +DH/2. 0 
Y B=YF-T-DH*DFLOAT (I-1) +DH/2.0 
EA=~YA*FY 
EB=--YB*FYF 
C=C+ (STRES (EA,SIGY,FYP,H,SIGU) ~STRES (EB, SIGY,FYP,H,SIGU) ) *DA 
S (I) =C/W/AH 
8 CONTINUE 
SG=0.0 ; 
DO 9 I=1,10 
U=GAM (S (I) ) 
UU=GAM (S (I+1) ) 
SG=SG+ (S (I) * (2. 0*U+UU) +S (I+1) *(2.0*UUtU) ) *DA/6.0*AH 
9 CONTINUE 
99 FORMAT ((5X,10F12.4)) 
RETURN 
END 


SUBROUTINE SNA (A,B,NN, NEQ) 
TPL LGM PRE Ane (A= 7 O-Z) 
DIMENSION A(NN,NN),B(NN) 
TF(NEQS EO: 1),\GO 7 TO) A 
NL=NEQ-1 
DO P3200 NN=1, NL 
IF (A(N,N).LE.0.0)GO TO 500 
N1=N+1 
DO 100 J=N1,NEQ 
100 A(N,J) =A(N,J) /A (N,N) 
DO 250 I=N1,NEOQ 
TFE(A(NJI)P.EQSO020)GO TO 250 
DO 200 J=1,NEQ 
NGL =A (i, J) A (LN) XA (Nd) 
200 Rip l Web Cl 73) 
250 B(I)=B (1I)—A (N,1) *B (N) 
300 B(N)=B(N) /A (N,N) 
M=NEQ 
B (MN) =B(M) ZA (M,M) 
DO 400 N=1,NL 
M1=M 
M=M-1 
pO 400 J=M1,NEQ 
400 B(M)=B(M)-B (J) *A (M, J) 
GO TO 600 
500 WRITE (6,1000) 
(CUMUAL, JE K IESE 
1000 FORMAT (' ZERO OR NEGATIVE ELEMENT ON NAIN DIAGONAL OF TRIANGURIZ 
1ED MATRIX') 
600 RETURN 
1 B(1) =B(1) ZA (1,1) 
GORTO R600 
END 
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SUBROUTINE MC 

IMPLICIT REAL*8 (A-H,O~2Z) 
COMMON/BL1/N,N1,AL,NM, NB (10) 
COMMON/TEA4/BC,TC,FC,EN,KT 
COMMON/TEA3/R2A,YRA, HR, RSNR, FYPR, DRS, RK 
COMNON/TEA2/SIGY,YS,H,RSN,RES,FYP,EPM,SIGU 
COMMON/TEA1/B,T, D,W,K1,K2, KF 
COMMON/BLS/SP (5) ,SCR(6,50,4) , PRO (4,100) 
NRA=0 

NRI=100 

NMI=NM+1 

ELM=3.0*10.**4 

IF (NM.EO.0)GO TO 20 

DOW 1 1=170N 1 

IF (I.NE. 1) JJ=NR (I-1) +1 

Dt Sf (| Chea W a FRc | 

IF (I. EQ.NM1) JJJ=N4+1 

IF (I.NE.NM1) JJJ=NR (1) 

JeJd0~Jd 

IF(J.GE.NRA) NRA=Jd 

IF (J. LE. NRI) NRI=J 

CONTINUE 

HH=AL/DFLOAT (N) 
SLA=HH*DFLOAT (NRA) /2.0 
SLI=HH*DFLOAT (NRI) /4.0 
IF(NM.EQ.0.0) SLA=AL 

IF (NN.EQ.0.0)SLI=AL/4.0 

DO 2 r= 175 

SP (I) =SLI+ (SLA~SL1I) *DFLOAT (I-1) /4.0 
RC=TCFEC 

AS=D¥*¥W+2.0*B¥*T 

Y 1= (AS+D/2. 0-AC¥TC/16. 0) / (AStAC/8. 0) 
FY=SIGY/ELM/D/2.0 

DOT sea 35 

po’ & y=1',5 

CALL NAPB(FY,BM,Y1,YD,SP(1) ) 

SCH (10s 1) =BM 

SCR (lL, Uy 2peary 

SCR (0,3) =¥1 

SCR Cl, J /4)=YD 

CONTINUE 

FY=1. 15*FY 

CONTINUE 

Y1= (AS*D/2. 0+RA* (D+DRS) ) / (AS+RA) 
FY=-SIGY/ELM/D/2.0 

po & =, £0 

CALI WARE CEY, BM YY by, YD 1.0) 

SCR (6,J,1)=BM 

SCR (6,7 2) =FY 

SCR (6, J, 3)= 11 

SCR(6,J,4) =¥D 

FY=1. 15*FY 

CONTINUE 

po?’ Gtr=S14,5 

WRITE (6,10) SP (T) 

pO" 7 1735 

WRITE (6,171) (SCR (I,J,K) ,K=1,4) 
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6 CONTINUE 
WRITE (6, 12) 
DOG 3=1,H0 
13 WRITE(6,11) (SCR (6,d,K) ,K=1,3) 
RETURN 
10 FORMAS (1H1,5X,'%*#*** POSITIVE MOMENT REGION *****', 8X, 
1@*#OSHEAR SPAN=',F7.2,." &*,//5%, (MOMENT! ,14X, 
2'*CURVATURE',11X,"*NUETRAL AXIS',8X,'2 N.A. DISTANCE!) 
12 FORNAT(1H1,5X,'***** NEGATIVE MOMENT REGION ****x"t, 
1//5X,* MOMENT", 14X,"CURVATURE', 11X, *‘NUETRAL AXIS") 
12 FORMAT (5X, 5(E15.7,5X)) 
END 


SUBROUTINE MMI (BM) 
IMPLICIY REAL*8 (A-H,O-Z) 
COMMONYBL1/N,N1,AL,NM,NR (10) 
CONNON/BL2/CIP,CIN, AKP, AKN, ATP, ATN, BYP, BYN 
COMMON /BL5S/SP (5) ,SCR(6,50,4) ,PRO(4,100) 
DIMFNSION BM (100),S (10) 
ELM=3.0*10.**4 
NM1=NM41 
CALL SSP(S,BM) 
N1=N+1 
DO 1 I=1,N1 
Js=1 
S PCE. (a9 2,843 
3 DSL=AL 
CI=CIP 
NR (NM1)=N1 
‘5 WES LESNRIS)) CGO TO 6 
EF (I.GE.NR (JS) ) JS=JS+1 
GO TO 5 
6 DO 4 IT=1,5 
D1=S (JS) -SP (IT) 
D1=DABS (D1) 
IF(D1-DSL) 20,20,4 
20 DST=p4 
JM=IL 
& CONTINUE 
1% K=1 
12 SM=BM (I) -SCR(JM,K,1) 
IF (EM (I) .LT.0.0) SM=-SM 
TCS, 7.8 
7 CCR SGT. 1) GO To 13 
BROM2).1) =SCR (3M, 1,2) *BM(1) /SCR(IM, 1, 1) 
PRO 3 AS) SSER (IM, 1,3) 
PRO (4,1) =SCR(JM,1,4) 
BROOCH eC PRO (2,1) -BN (1) 
Go TO 1 
8 K=K+1 
IF(K.EQ.40) STOP 100 
GOlleo.G12 
13 IF(Bi (1) .LT.0.0) SM=-SM 
po 14 Ir=2,4 
16 (PROG TAT) =SCR (JN, K,1L) + (SCR (OM, K, 11) -SCR (OM, (K=1) , 11) ) ¥54/ 
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’ 


T(SCR(JM,K, 1) -SCR (JM, ( 
Crs 


PRO(1,I)= CI *PRO 
éotTo «4 

2 JM=6 
Creer 
Gortoe11 

1 CONTINUE 

' 98 FORMAT (5X, 515) 

99 FORMAT (5X, 8E15.5) 

END 


SUBROUTINE DERM(Q,P,BM,SH,R,DEF,SR,SD) 
IMPLICIT REAL*8 (A-H,O-Z) 
COMMON/BL1/N,N1,AL,NM,NR (10) 
COMMON/BL2/CIP,CIN, AKP, AKN, ATP, ATN, BYP, BYN 
COMMON/BL5/SP (5) ,SCR(6,50,4) , PRO (4,100) 
DIMENSION X(10) ,XX(10).,M(4) ,Q(100) ,2 (100) , BM(100) , EI (100) ,G (100) 
1,R (100) ,EM (100) ,DEF (100) ,SH(100) ,SR(100) ,SD(100) ,X3(10) 
H=AL/DFLOAT (N) 
CX=1.0 
CP=1.000001 
A3=30.0 
IF(NM.EQ-.0)GO TO 41 
DO 1 I=1,NM 

£ Ve (2-070 

44 ELM=3.0*10.0**4 
DO 40 I=1,N1 

4G EM(ZI)=0.0 
M(1)=1 
M(2)=1 
M (3) =1 
M (4) =0 
MR=0 
POC AM TEAC ANA 
BE (i) CIR 

14 CONTINUE 
IF(NM.EQ.0)GO TO 42 

7 CALUPREACT (1,0,P,EM, BM, B1,G,X) 
pO 3 I=1,NM 
P (NR (I))=XX (I) + (X (I) -XX (I) ) /A3 
IF(MR.LE. 1) P (NR (I) ) =X (I) 
IF (MR. EQ. 2) P (WR (I) ) =X (I) *0. 98 

3 CONTINUE 
CA=01.0 
€B=0...0 
pO 8 I=1,NN 
CA=CA+ (X (I) ~XX (I) ) **2 
CB=CB+X(1) **2 

8 CONTINUE 
IF (MR.LT.10) AA=0.005 
TF(MR.GE.10) AA=0.01 
IF (MR.GT.30)STOP 10 
CC=DSQRT(CA/CB) 
ir cch it Ceyco, Toren 
DO 62 I=1,NM 
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62 P(NR(1)) =X3 (1 
coo iy 


61 IF(CC.LT.CP) CP=CC+0.00001 
WRITE (6,51) CC, (X(I) ,1=1, NM) 
51 FORMAT (5X,'**** ERROR ESTIMATION IS 
T10F7. 2) 
S39 I P(CCeLELAA) GO TO 22 
42 CALL FDM(M,Q,P,EM,BM,EL,G,SH,R,DEF) 
DO 30 I=1,N1 
30 BM (I) =BM (I) -EM (1) 
DO 4 I=1,N1 
IF (BM(I))5,5,6 
5 EI (I)=CIN 
GO TO 4 
60FT (1) =CIP 
4 CONTINUE 
IF(NM.EQ.0)GO TO 45 
_ IF(MR.LE.1) GO TO 12 
“5 CALL MMI (BM) 
DO ou T=1,81 
44 EM (I) =PRO (1,1) 
IF(NM.EQ.0)GO TO 43 
DO 9 I=1,NM 
So Cl) =xe (0) 
9 XX (I) =P (NR (I)) 
12 MR=MR+1 
60 DO 15 I=1,NM 
15 P(NRE(I)) =0.0 
GOTO. 7 
22,0034 Te1, NM 
SIP (NE (2) (2) 
43 CALL FDM(M,Q,P,EM,BM,EL,G,SH,R,DEF) 
DOss28 1,01 
32 BN (I) =BM (i) -EM (I) 
- CALL DS (SH, PRO,SD) 
DO) 23) ie 1,N1 
23 SD(I) =DEF(I)+SD (I) 
DO 21 1=2,8 
21 SR(L)=(SD(I+1)-SD(I-1)) /H/2.0 
SR(1) =(SD(2)-SD(1)) /H 
SR(N1)=(SD(N1)-SD(N)) /H 
100 FORMAT (5X, 5E15.7) . 
RETURN 
99 FORMAT ((3X, 10E12.3)) 
END 
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